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ABSTRACT 

The  Transverse  Ranges  province  consists  of  a  set  of  east-trending  mountains  and  geologic  structures  that  are  distinct 
from  the  general  northwest-southeast  groin  of  the  other  provinces  of  California.  The  Coast  Ranges  province  to  the  north 
and  the  Peninsular  Ranges  province  to  the  south  are  characterized  by  strike— slip  movement  on  reasonably  straight  faults 
striking  northwest.  In  the  Transverse  Ranges,  thrust  faulting  is  common.  Most  of  the  strike-slip  faults  are  curved  and  strike 
approximately  east.  The  persistent  northwesterly  trend  of  the  Son  Andreas  fault  changes  to  on  east  trend  in  the  Transverse 
Ranges  forming  a  bend  in  the  fault. 

The  proximity  of  the  bend  in  the  Son  Andreas  fault  to  the  Transverse  Ranges  suggests  that  there  may  be  a  direct 
relationship  between  the  San  Andreas  fault  and  the  structures  in  the  Transverse  Ranges.  In  this  work,  the  bend  in  the  San 
Andreas  fault  is  approximated  by  a  series  of  straight  dislocation  surfaces,  and  the  displacement  and  stress  fields  produced 
at  an  observation  point  by  each  segment  are  summed.  Two  models  with  "locked"  segments  on  the  Son  Andreas  fault  15 
and  30  km  deep  are  studied. 

The  calculated  vertical  displacement  fields  from  the  two  models  resemble  the  general  topographic  features  of  the 
Transverse  Ranges,  suggesting  that  the  present  topography  may  hove  been  produced  by  the  bend  in  the  San  Andreas 
fault.  It  is  not  possible  to  select  between  the  two  models  of  the  San  Andreas  fault  based  on  the  topography  alone. 

Several  conclusions  result  from  the  interpretation  of  the  stress  fields  calculated  from  the  two  models. 

First,  it  does  appear  to  be  possible  to  select  between  the  two  models  based  on  the  calculated  stress  fields.  Model  2, 
the  model  with  the  locked  zone  30  km  deep,  seems  to  explain  the  observed  fault  patterns  and  strain  data  better  than 
model  1,  the  model  with  the  locked  zone  15  km  deep.  If  this  model  is  correct,  the  sliding  port  of  the  fault  is  deeper  in 
the  Transverse  Ranges  than  in  central  Colifornia. 

Second,  the  principal  compressive  stresses  calculated  from  the  model  have  a  north-south  orientation  throughout  most 
of  the  study  area.  This  agrees  with  the  conclusions  drawn  by  other  workers. 

Third,  the  areas  of  potentiol  failure  determined  from  the  model  include  mony  of  the  larger  earthquakes  in  the  Transverse 
Ranges;  and  the  larger  events  that  have  not  occurred  in,  or  adjacent  to,  the  area  of  potential  failure  can  be  related  to 
movement  on  faults  in  the  Peninsular  Ranges. 

Finally,  the  calculated  stress  system  could  explain  the  pattern  of  the  faults  in  the  Transverse  Ranges  having  Holocene 
or  Quaternary  movement.  This  does  not  necessarily  imply  that  the  bend  in  the  San  Andreas  fault  is  responsible  for  the 
origin  of  the  Transverse  Ranges.  It  does,  however,  mean  that,  given  a  complex  set  of  faults  of  many  different  orientations, 
the  stress  system  produced  by  the  bend  in  the  San  Andreas  fault  will  cause  activity  mainly  on  those  faults  that  are  favorably 
oriented  to  the  stress  system  at  the  time  of  observation.  For  the  some  time  interval,  other  less  favorably  oriented  faults 
will  show  a  much  lower  level  of  activity,  even  though  these  faults  may  hove  been  much  more  active  at  some  earlier  time. 
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INTRODUCTION 

It  is  generally  accepted  among  earth  scientists  that  the 
surface  of  the  earth  can  be  divided  into  several  plates  that 
are  moving  relative  to  one  another  (Morgan,  1968;  Le 
Pichon,  1968;  Atwater  and  Molnar,  1973).  Each  plate  has 
a  source  (generally  a  mid-ocean  rise)  and  a  sink  (either 
a  trench  or  a  mountain  range).  The  plates  can  move  past 
one  another  along  transform  faults. 

Atwater  (1970)  and  Anderson  (1970)  suggested  that 
many  of  the  geologic  features  of  California  could  be  ex- 
plained in  terms  of  the  interaction  between  the  Pacific  and 
North  American  plates.  They  suggested  that  the  San  An- 
dreas fault  is  behaving  as  a  transform  fault  to  accommo- 
date most  or  all  of  the  movement  between  the  two  plates. 
However,  Hill  (1971,  1974)  has  strongly  challenged  this 
interpretation. 

The  transform  faults  mapped  in  the  ocean  floors  are 
generally  very  linear  features  (Morgan,  1968;  Le  Pichon, 
1968)  while  the  San  Andreas  fault  trace  contains  several 
large  bends.  The  bend  in  the  San  Andreas  fault  near  Tejon 
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Pass  (figure  1 )  is  of  particular  interest  to  California  geolo- 
gists because  of  its  proximity  to  the  Transverse  Ranges 
geomorphic  province. 

The  Transverse  Ranges  province  consists  of  a  distinct 
set  of  east-trending  mountains  and  geologic  structures 
that  cut  across  the  general  northwest-trending  grain  of 
the  other  provinces  of  California  (Bailey  and  Jahns,  1954, 
figure  1 ) .  The  faulting  of  the  Coast  Ranges  province  to  the 
north  (Page,  1966)  and  the  Peninsular  Ranges  province 
to  the  south  (Jahns,  1954)  is  primarily  characterized  by 
strike-slip  movement  on  reasonably  straight  faults  strik- 
ing northwest  with  only  minor  dijj-slip  activity.  However, 
in  the  Transverse  Ranges,  there  is  a  great  deal  of  thrust- 
faulting  (Bailey  and  Jahns,  1954;  Dibblee,  1968;  Crowell, 
1968).  Most  of  the  strike-slip  faults  are  curved  and  strike 
approximately  east.  The  persistent  northwesterly  trend  of 
the  San  Andreas  fault  in  northern  California  changes  to 
an  east  strike  in  the  Transverse  Ranges  forming  a  bend  in 
the  fault. 

The  origin  of  the  structures  in  the  Transverse  Ranges 
is  still  subject  to  much  debate.  There  is  some  indication 
that  the  east-trending  structures  may  have  developed  by 
late  Cretaceous  or  early  Tertiary  time  (Baird  and  others, 
1974);  but  the  Transverse  Ranges  trends  were  clearly  de- 
veloped by  Miocene  time  and  the  present  structures  have 
developed  since  Miocene  time  (Jahns,  1973). 

Benioff  (1955),  Jahns  (1973),  and  many  others  have 
suggested  that  the  bend  in  the  San  Andreas  fault  might  be 
responsible  for  the  post-Miocene  compressive  structures 
of  the  Transverse  Ranges.  Intuitively,  this  is  reasonable, 
for  one  would  expect  compressive  structures  on  the  con- 
cave sides  of  the  bends  in  the  San  Andreas  fault.  Addition- 
ally, there  seems  to  be  reasonably  good  evidence  that  part 
of  the  San  Andreas  fault  between  Tejon  Pass  and  Cajon 
Pass  is  no  older  than  Miocene  (Jahns,  1973;  Baird  and 
others,  1974;  Berry,  1974).  Thus,  the  compressive  struc- 
tures of  the  Transverse  Ranges  and  the  severe  bend  in  the 
San  Andreas  seem  to  be  about  the  same  age. 
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DEFORMATION.  STRESS  AND  FAULTING  TRANSVERSE  RANGES 


The  current  tectonic  activity  of  the  Transverse  Ranges 
is  not  clear.  There  is  a  history  of  moderate  seismicity 
(Hileman  and  others,  1974),  and  the  1971  San  Fernando 
earthquake  serves  as  a  reminder  that  the  province  is  still 
active.  Most  of  the  east-west  trending  faults  seem  to  have 
been  active  in  the  last  1  to  2  million  years,  and  some  of 
these  faults  have  clear  evidence  of  Holocene  movement 
(Jennings,  1973). 

The  proximity  of  the  bend  in  the  San  Andreas  fault  to 
the  Transverse  Ranges  suggests  that  there  may  be  a  direct 
relationship  between  them.  In  this  paper,  results  from 
elastic  dislocation  theory  (Chinnery,  1961,  1963)  will  be 
used  to  develop  a  mathematical  model  of  the  San  Andreas 
fault  in  southern  California.  The  deformation  and  stress 
fields  calculated  from  the  model  will  be  compared  to  the 
topography  and  the  mapped  geology  to  determine  any 
relationships  between  the  San  Andreas  and  the  Transverse 
Ranges  area.  The  boundaries  of  the  study  area  are  larger 
than  those  of  the  Transverse  Ranges  geomorphic  prov- 
ince. 

STRUCTURAL  CHARACTERISTICS 
San  Andreas  Fault 

The  San  Andreas  fault  can  be  traced  continuously  for 
520  miles  from  Point  Arena  in  northern  California  to  San 
Gorgonio  Pass  in  southern  California  (Crowell,  1962,  fig- 
ures 3  and  5).  Allen  (1957)  was  unable  to  trace  the  fault 
through  San  Gorgonio  Pass;  rather,  he  shows  it  dying  out 
before  it  reaches  the  Banning  fault.  He  suggests  that  the 
San  Jacinto  fault  is  now  the  active  trace  of  the  San  An- 
dreas fault.  However,  Dibblee  (1968)  shows  the  San  An- 
dreas fault  branching  north  of  San  Gorgonio  Pass  and 
continuing  southward.  If  the  continuity  of  the  San  An- 
dreas fault  through  the  San  Gorgonio  Pass  area  is  accept- 
ed, the  fault  may  be  traced  nearly  100  miles  farther  south 
to  the  eastern  side  of  the  Salton  Sea.  There,  it  is  lost  under 
alluvium.  If  present  seismicity  is  accepted  as  a  reasonably 
good  indication  of  the  location  of  the  fault  trace,  then  the 
broad  line  of  earthquake  epicenters  running  from  the  Sal- 
ton  Sea  south  to  the  Gulf  of  California  (figure  3)  suggests 
that  the  fault  zone  continues  on  into  the  Gulf  (Richter, 
1958,  figure  28-1). 

In  northern  California,  the  fault  goes  out  to  sea  at  Point 
Arena.  However,  ground  breakage  occurred  at  Shelter 
Cove  after  the  1906  earthquake,  and  Curry  and  Nason 
(1967)  show  the  fault  curving  north  parallel  to  the  coast 
up  to  Shelter  Cove  rather  than  continuing  northwest  to 
Point  Arena  to  explain  both  the  ground  breakage  at  Shel- 
ter Cove  and  the  seismic  profiles  northwest  of  Point  Arena 
that  show  no  evidence  of  the  northwest  extension.  On  the 
basis  of  earthquake  epicenters,  magnetic  anomalies,  and 
topography,  Oakeshott  ( 1966)  and  Nason  (1968)  suggest 
that  the  fault  zone  bends  at  Cape  Mendocino  to  strike 
nearly  east-west  along  the  Mendocino  Escarpment.  On 
the  basis  of  very  accurate  epicenter  determinations  and 
first-motion  studies.  Bolt  and  others  (1968)  suggest  that 
the  fault  divides  into  two  branches  near  Cape  Mendocino. 


One  branch  swings  west  to  join  the  Mendocino  Escarp- 
ment, and  the  other  branch  continues  on  a  northwesterly 
strike  from  Cape  Mendocino.  If  the  extension  of  the  fault 
northward  as  far  as  the  Mendocino  Escarpment  is  accept- 
ed, the  San  Andreas  fault  is  well  over  100  miles  long, 
running  from  Cape  Mendocino  in  the  north  to  the  head 
of  the  Gulf  of  California  in  the  south. 

North  and  south  of  the  Transverse  Ranges,  the  strike  of 
the  San  Andreas  fault  is  very  nearly  northwest-southeast, 
but  in  the  Transverse  Ranges  the  strike  becomes  west  of 
northwest  (figure  2).  The  dip  of  the  fault  in  the  study  area 
is  approximately  vertical,  and  the  sense  of  motion  is  right 
lateral.  The  amount  of  displacement  on  the  fault  in  the 
study  area  appears  to  be  160-175  miles  since  the  Oligo- 
cene  (Crowell,  1962).  Sag  ponds  and  offset  stream  chan- 
nels indicate  that  the  San  Andreas  fault  has  been  active  in 
the  Holocene  in  the  study  area.  However,  recent  surveys 
of  small  geodetic  figures  crossing  the  fault  from  1964- 
1972  show  no  movement  on  the  San  Andreas  fault  (Ap- 
pendix B).  Allen  (1968)  suggests  that  fault  activity  on 
this  part  of  the  San  Andreas  fault  is  evidenced  by  large 
infrequent  earthquakes  rather  than  by  small  earthquakes 
and  fault  creep. 

The  level  of  seismic  activity  is  high  along  all  of  the  San 
Andreas  fault  system  from  Hollister  to  the  Mexican  bor- 
der, except  on  the  fault  segment  from  the  Carrizo  Plain  to 
the  vicinity  of  Cajon  Pass  ("Seismicity"  in  this  report; 
Bolt  and  others,  1968;  Bolt  and  Miller,  1971;  Allen  and 
others,  1965;  Hileman  and  others,  1973).  A  microearth- 
quake  survey  of  the  San  Andreas  fault  system  also  sup- 
ports this  picture  of  the  seismic  activity  (Brune  and  Allen, 
1967).  However,  there  is  evidence  of  historic  movement 
on  this  segment  of  the  fault.  The  1857  Fort  Tejon  earth- 
quake apparently  broke  the  ground  surface  from  the  vicin- 
ity of  Cholame  to  somewhere  between  Cajon  Pass  and  San 
Gorgonio  Pass  (Wood,  1955;  Allen,  1968).  These  obser- 
vations have  led  to  the  suggestion  that  this  section  of  the 
fault  is  "locked"  and  is  a  likely  spot  for  the  next  great 
earthquake  in  California  (Allen,  1968). 

San  Gabriel  Fault 

The  San  Gabriel  fault  lies  to  the  west  of  the  San  An- 
dreas fault  in  the  Transverse  Ranges.  The  northern  part 
of  the  fault  trends  approximately  northwest,  while  the 
southern  part  of  the  fault  has  a  more  easterly  trend.  Cro- 
well (1954)  has  presented  evidence  of  about  20  miles  of 
right-lateral  strike-slip  displacement  since  the  early  Mio- 
cene on  the  San  Gabriel  fault.  There  is  some  evidence 
(Crowell,  1968)  that  the  southern  part  of  the  San  Gabriel 
fault  has  been  active  in  the  Holocene,  but  there  is  no 
evidence  of  Holocene  activity  on  the  northern  section. 
However,  Jennings  (1973)  indicates  Quaternary  displace- 
ment on  the  northern  part  of  the  fault.  Jahns  (1973) 
suggests  that,  in  the  Miocene,  the  San  Gabriel  fault  was 
the  main  trace  of  the  San  Andreas  fault  system.  Jahns' 
reconstruction  indicates  that,  by  the  middle  Pliocene,  the 
San  Andreas  bend  had  formed  and  the  San  Gabriel  fault 
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was  no  longer  accounting  for  a  major  amount  of  the  mo- 
tion on  the  San  Andreas  system. 

The  San  Andreas  and  the  San  Gabriel  faults  are  the 
major  right-lateral  strike-slip  faults  in  the  Transverse 
Ranges.  Since  Jahns'  reconstruction  suggests  that  the 
bend  of  the  San  Andreas  fault  did  not  form  until  middle 
Pliocene  time,  in  this  paper  I  will  concern  myself  only  with 
the  geologic  history  since  middle  Pliocene  time. 

Garlock  Fault 

The  Garlock  fault  is  a  northeast-striking  left-lateral 
strike-slip  fault.  The  fault  is  approximately  260  km  long, 
and  evidence  can  be  found  for  38  to  64  km  of  left-lateral 
displacement  (Davis  and  Burchfiel,  1973).  There  is  also 
some  evidence  from  geodetic  data  (Greensfelder  and  Ben- 
nett, 1973;  Appendix  B)  and  from  offset  stream  channels 
and  deposits  that  the  Garlock  fault  is  currently  active. 
Davis  and  Burchfiel  ( 1973)  presented  a  summary  of  what 
is  known  about  the  Garlock  fault  together  with  an  inter- 
esting hypothesis  about  its  origin.  They  suggest  that  the 
Garlock  fault  may  be  an  intracontinental  transform  fault 
rather  than  a  major  transcurrent  fault,  and  their  argu- 
ments are  fairly  convincing. 

White  Wolf  Fault 

The  White  Wolf  fault  is  a  left-lateral  strike-slip  fault 
parallel  to,  and  northwest  of,  the  Garlock  fault.  It  is  a 
high-angle  reverse  fault,  dipping  to  the  southeast  with  the 
southern  block  moving  up  relative  to  the  northern  block. 
The  White  Wolf  fault  is  active  and  was  the  locus  of  the 
1952  Arvin-Tehachapi  earthquakes  (Oakeshott,  1955). 

On  July  21,  1952,  at  4:52  AM  PDT,  a  7.7  magnitude 
earthquake  occurred  on  the  White  Wolf  fault  and  was 
accompanied  by  17  miles  of  surface  faulting.  Mapping  of 
these  movements  indicated  up  to  2  feet  of  left-lateral 
strike-slip  movement  and  2  to  3  feet  of  vertical  motion. 
The  instrumental  fault  plane  solution  showed  that  the 
strike  of  the  fault  plane  was  about  north  50°  east  and  the 
dip  was  60°  to  66°  to  the  southeast  (Gutenberg,  1955).  The 
fault  plane  solution  also  indicated  that  the  south  side  of 
the  fault  moved  up  relative  to  the  north  side  and  that  there 
was  left-lateral  strike-slip  motion;  thus,  the  fault  plane 
solution  agreed  with  observed  surface  faulting.  Triangula- 
tion  and  leveling  surveys  after  the  earthquake  by  the  U.S. 
Coast  and  Geodetic  Survey  showed  that  there  was  about 
2  feet  of  left-lateral  horizontal  offset  and  3  to  4  feet  of 
vertical  offset,  with  the  south  side  of  the  fault  moving  up 
relative  to  the  north  side  (Oakeshott,  1955). 

BeniofT  (1955,  p.  204)  suggested  that  "...the  distortion 
of  that  portion  of  the  eastern  block  of  the  San  Andreas 
fault  moving  southward  along  the  curved  restraint  of  the 
great  bend  must  be  principally  in  the  form  of  a  compres- 
sion oriented  approximately  north-south.  The  White 
Wolf  fault  is  thus  a  mechanism  for  relief  of  this  localized 
stress." 


Big  Pine  Fault 

The  Big  Pine  fault  strikes  approximately  east-west  and 
extends  for  about  50  miles  to  the  west  of  the  San  Andreas 
fault.  Its  topographic  expressions  are  similar  to  those  of 
the  San  Andreas,  the  Garlock,  and  other  large  strike-slip 
faults.  There  may  be  as  much  as  14  miles  of  left-lateral 
offset  (Hill,  1955;  Carman,  1964);  and  in  the  Lockwood 
Valley  area,  there  is  some  evidence  of  dip-slip  motion, 
with  the  north  side  of  the  fault  moving  up  relative  to  the 
south  side  on  an  approximately  vertical  fault  (Carman, 
1964).  There  is  evidence  that  the  Big  Pine  fault  was  active 
as  recently  as  the  Pleistocene;  but  there  is  no  conclusive 
evidence  that  it  has  been  active  in  the  Holocene  (Carman, 
1964),  although  a  large  earthquake  that  occurred  in  1852 
may  have  had  its  epicenter  on  the  Big  Pine  fault  (Townley 
and  Allen,  1939).  Geodimeter  surveys  from  about  1962  to 
the  present  do  not  indicate  horizontal  movement  (Hof- 
mann,  1968;  California  Division  of  Mines  and  Geology, 
unpublished  data,  1973). 

Santa  Ynez  Fault 

The  Santa  Ynez  fault  is  an  east-west  trending  fault 
lying  to  the  west  of  the  San  Andreas  fault.  It  has  been 
described  in  some  detail  by  Dibblee  (1950,  1966).  The 
fault  dip  is  vertical  along  the  western  section  of  the  fault, 
approximately  50°  to  the  south  along  the  central  section 
of  the  fault,  and  very  steeply  dipping  to  the  north  along 
the  eastern  section  of  the  fault.  Dibblee  (1950,  1966) 
documented  several  miles  of  vertical  offset  on  the  Santa 
Ynez  fault,  the  south  side  of  the  fault  being  elevated  rela- 
tive to  the  north  side,  and  also  found  some  evidence  of 
horizontal  offset,  with  the  offset  being  in  a  left-lateral 
sense.  The  fault  was  certainly  active  in  the  Pleistocene  and 
possibly  in  the  Holocene  (Dibblee,  1966,  p.  76).  A  recent 
analysis  of  a  first-order  level  survey  suggests  that  there 
was  no  vertical  movement  on  the  fault  from  1920  to  1971 
(Buchanan  and  others,  1973).  Schroeter  (1972)  presents 
some  evidence  that  the  central  part  of  the  Santa  Ynez  fault 
system  has  been  active  in  the  Holocene. 

Thrust  Faults 

The  study  area  contains  many  thrust  faults  (Hill,  1955, 
figure  1;  Bailey  and  Jahns,  1954,  plate  4;  Bailey,  1954; 
Merrill,  1954;  Crowell,  1968). 

The  Pleito  fault,  north  of  the  San  Andreas  fault,  seems 
to  have  been  active  in  the  Holocene  (Crowell,  1968).  It 
dips  shallowly  to  the  south,  and  the  south  side  has  moved 
up  relative  to  the  north  side. 

The  Frazier  Mountain  and  Pine  Mountain  thrust  faults, 
south  of  the  San  Andreas  fault,  were  active  in  the  Pleisto- 
cene (Crowell,  1968);  but  there  does  not  appear  to  be  any 
reliable  evidence  of  activity  in  the  Holocene.  As  can  be 
seen  on  figure  2,  motion  on  the  Frazier  and  the  Pine 
Mountain  faults  is  to  the  south. 

The  San  Cayetano,  the  Oak  Ridge,  the  Santa  Susana, 
and  the  Sierra  Madre  faults  may  have  been  active  in  the 


CALIFORNIA  DIVISION  OF  MINES  AND  GEOLOGY 


BULL.  203 


HolcKene  (Jennings.  1973).  Again,  figure  2  shows  that 
motion  on  the  San  Cayetano,  the  Santa  Susana.  and  the 
Sierra  Madrc  faults  is  to  the  south,  while  motion  on  the 
Oak  Ridge  fault  is  to  the  north. 

Finally,  the  zone  of  faulting  defined  by  the  Santa  Mon- 
ica, the  Raymond  Hill,  and  the  Sierra  Madre  faults  is 
characterized  by  a  left-lateral  displacement  and  seems  to 
have  been  active  in  the  Holocene  (Yerkes  and  others, 
1965). 

SEISMICITY 

The  seismicity  of  southern  California  has  been  dis- 
cussed in  detail  by  Richter  ( 1958)  and  Allen  and  others 
( 1965).  and  Hileman  and  others  (1973)  published  a  cata- 
log of  instrumentally  located  earthquakes  in  southern  Cal- 
ifornia from  1932  to  1972.  Figure  3  shows  that  the  San 
Andreas  fault  north  of  the  junction  with  the  San  Jacinto 
fault  has  been  relatively  inactive,  while  the  San  Andreas 
has  been  moderately  active  between  the  San  Jacinto  junc- 
tion and  the  Salton  Sea.  The  major  seismic  activity  in 
southern  California  apparently  was  associated  with  the 
San  Jacinto,  the  Newport-Inglewood.  and  the  White  Wolf 
faults.  The  small  area  near  the  San  Gabriel  fault  is  the 
aftershock  zone  of  the  1971  San  Fernando  earthquake. 

Comparing  figures  4  and  5.  it  is  clear  that  not  all  large 
historic  events  have  documented  surface  faulting.  Howev- 
er, there  is  some  evidence  of  surface  faulting  associated 
with  the  1852  earthquake  on  the  Big  Pine  fault  (Townley 
and  Allen.  1939;  Carman.  1964).  There  is  also  some  indi- 
cation of  surface  breakage  associated  with  either  the  1918 
or  the  1923  earthquakes  on  the  San  Jacinto  fault  (Donald 
L.  Fife.  California  Division  of  Mines  and  Geology,  person- 
al communication.  1974). 

The  1857  Fort  Tejon  earthquake  on  the  San  Andreas 
fault  was  one  of  the  three  largest  earthquakes  in  Califor- 
nia's recorded  history  (Wood,  1955),  and  the  length  of  the 
fault  break  is  well  documented.  There  is  reasonably  good 
e\  idence  that  surface  breakage  extended  to  the  vicinity  of 
San  Bernardino  (Townley  and  Allen,  1939;  Wood,  1955; 
California  Division  of  Mines  and  Geology,  unpublished 
data),  and  there  is  tenuous  evidence  that  the  faulting  ex- 
tended to  the  south  of  San  Bernardino  (Wood,  1955).  The 
northern  end  of  the  fault  break  was  probably  in  the  vicin- 
ity of  Cholame  (Wood,  1955).  The  descriptions  of  the 
surface  breaks  clearly  indicate  large  right-lateral  offset, 
but  the  amount  of  offset  was  not  measured  anywhere. 
However,  indirect  evidence  from  offset  stream  channels  in 
the  Carrizo  Plain  suggests  that  the  offset  may  have  been 
as  much  as  10  meters  (Wallace,  1968). 

Relatively  few  of  the  earthquakes  on  figure  5  are  located 
north  of  the  trend  defined  by  the  Malibu  Santa  Monica- 
Raymond-Sierra  Madre-Cucamonga  faults.  This  trend  is 
the  traditional  southern  boundary  of  the  Transverse 
Ranges  province,  and  the  faults  south  of  the  trend  have 
the  same  northwest  orientation  as  the  San  Andreas  sys- 


tem. In  fact,  an  argument  can  be  made  that  faults  such  as 
the  Newport-Inglewood  and  the  Elsinore  are  in  some  way 
ancestral  to  the  San  Andreas  fault  (Berry,  1973). 

The  magnitude  7.2  Arvin-Tehachapi  earthquake  in 
1952  was  the  largest  earthquake  in  the  study  area  since 
1857.  Its  effects  are  described  in  "Structural  Characteris- 
tics". Earthquakes  in  1952  and  1953  in  the  vicinity  of  the 
White  Wolf  and  the  Pleito  faults  may  have  been  after- 
shocks of  the  1952  earthquake,  but  the  occurrence  of  two 
earthquakes  in  that  area  before  1952  suggests  that  these 
post- 1952  earthquakes  might  not  have  been  aftershocks. 

The  earthquake  on  April  4,  1893.  between  the  San  Ga- 
briel and  the  San  Andreas  faults  could  have  occurred  on 
a  western  extension  of  the  fault  system  that  produced  the 
1971  San  Fernando  earthquake.  Newspaper  reports  indi- 
cate that  damage  was  greater  to  the  west  and  northwest 
of  San  Fernando  than  in  San  Fernando  (Townley  and 
Allen,  1939).  There  is  also  evidence  that  indicates  the 
weak  possibility  of  surface  faulting. 

The  series  of  earthquakes  at  the  western  end  of  the 
Malibu  fault  plus  a  magnitude  5.5  earthquake  in  1973 
(Ellsworth  and  others,  1973)  in  the  same  area  suggest  that 
the  Malibu  fault  may  be  still  active. 

Historically,  the  Santa  Barbara  Channel  has  been  the 
scene  of  minor  to  moderate  seismic  activity  (Hamilton 
and  others.  1969;  Sylvester  and  others.  1970;  Lee  and 
Vedder.  1973);  figure  3  shows  the  moderate  earthquakes. 
Sylvester  and  others  (1970)  discussed  an  earthquake 
swarm  that  occurred  in  the  Santa  Barbara  Channel  in 
1968.  Movement  on  a  northwest-striking  structure  is  the 
preferred  solution  according  to  Sylvester  and  others 
(1970).  although  the  major  structural  features  of  the 
Santa  Barbara  Channel  strike  approximately  east-west.  In 
contrast.  Lee  and  Vedder  (1973)  produced  a  composite 
fault  plane  solution  for  earthquakes  on  or  near  one  of  the 
faults  in  the  Channel,  and  they  found  left-lateral  move- 
ment on  an  east-west  plane. 

Whitcomb  and  others  ( 1973)  discuss  the  relationship  of 
the  1971  earthquake  to  the  regional  seismicity  pattern. 
They  suggest  that  the  San  Fernando  earihquakc  may  have 
occurred  on  a  zone  they  call  the  Dume  I'alnulale  zone 
and  that  this  zone  shows  up  in  seismicity  trends  from  I960 
to  1970  (Whitcomb  and  others.  1973.  figure  22).  Howev- 
er, study  of  their  epicenter  plot  (figure  22A)  indicates  that 
this  zone  is  only  vaguely  developed  in  the  data,  if  at  all. 
Their  epicenter  plot  suggests  rather  that  the  Transverse 
Ranges  -  Los  Angeles  Basin  area,  is  an  area  of  diffuse 
minor  seismicity  Maps  published  by  Hileman  and  others 
(1973)  also  support  this  interpretation 

GEODETIC  CHARACTERISTICS 

Since  I  he  1906  earihquakc  in  San  Francisco,  much  ef- 
fort has  gone  into  the  study  of  crustal  movement  using 
geodetic  techniques.  Greensfelder  (1972)  listed  the  loca- 
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Figure  3.        Epicenters  of  earthquakes  in  southern  California  wilf»  magnitude  4  and  greater,  1932-72   {from  Hileman  ond  others,  1973,  p. 
64). 
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Figure  6.        Triangulation  neli  for  crustal  movement  studies   (modified  from  Greensfelder,  1972.  p.  5). 
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Figure  7.        Geodiineter  fault-monitoring  program   (from  Greensfelder,  1972,  p.  11). 
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tion  and  dates  of  these  geodetic  studies.  Many  analyses  of 
the  data  have  been  made  (sec,  e.g.,  Burford,  1*^66,  and 
Savage  and  Burford,  1973,  for  partial  hsts  of  references). 

Triangulation  Surveys 
Taft-Mojave  Net 

In  1959-1960,  a  large  triangulation  network  was  estab- 
lished in  the  study  area  (figure  6)  by  the  U.S.  Coast  and 
Geodetic  Survey  and  the  California  Department  of  Water 
Resources  (Miller  and  others,  1969a).  Many  of  the  line 
lengths  in  the  net  were  measured  using  a  Model  2A  Geo- 
dimeter.  In  1967,  the  net  was  re-triangulated,  and  again 
selected  lines  were  measured  with  the  Geodimeter.  The 
following  results  are  quoted  from  the  report  (Miller  and 
others,  1969a,  page  ii): 

"The  resultoni  position  vectors  and  their  corresponding  95%  confidence 
error  ellipses  show  thot  the  chonges  in  position  ore  within  the  limits  of 
observotion.  However,  significant  patterns  exist.  Strain  is  more  evident  than 
fault  slippage.  There  is  an  area  of  expansion  east  of  the  Gorlock  fault  and 
north  of  the  San  Gabriel  fault.  The  remainder  of  the  net  generally  contracts 
or  compresses.  Along  the  White  Wolf  fault,  the  Geodimeter  results  show 
there  is  left-lateral  movement  of  about  9  cm  for  the  interval  1959-1960  to 
1967.  On  the  Son  Andreas  foult,  west  of  the  junction  with  the  Gorlock  ond 
San  Gabriel  faults,  the  Geodimeter  results  indicate  right-lateral  movement 
of  about  9  cm  for  the  7-yeQr  period  Along  the  Gorlock  fault  ond  the  San 
Andreas  fault  east  of  the  junction,  the  results  indicate  stroin  rather  than 
slippage." 


The  establishment  of  this  net  represents  a  major  step 
toward  the  study  of  crustal  deformation  in  the  study  area. 
However,  the  results  of  the  resurvey  cannot  be  accepted 
without  reservation.  First,  the  overall  accuracy  of  the  sur- 
vey was  about  4  parts  per  million  (Miller  and  others, 
1969a,  page  ii).  Thus,  the  results  obtained  are  well  within 
survey  error,  although  the  consistency  of  the  results  does 
give  them  greater  credibility.  Second,  and  perhaps  more 
important,  if  the  hypothesis  of  a  time-varying  strain  field 
advanced  by  Greensfelder  and  Bennett  (1973)  is  valid,  it 
may  be  impossible  to  separate  the  effects  of  the  time- 
varying  strains  and  the  tectonic  strains.  Certainly  more 
resurveys  of  the  net  will  be  necessary  to  substantiate  the 
results. 

San  Fernando-Bakersfield  Arc 

This  triangulation  arc  (figure  6)  was  established  in 
1932  and  resurveyed  in  1952  and  1962  (Meade  and  Small, 
1966).  Some  of  the  stations  in  the  arc  are  common  to  the 
Taft-Mojave  net  (Miller  and  others,  1969b),  and  an  ex- 
tensive analysis  of  the  1962  survey  is  in  essential  agree- 
ment with  the  results  from  the  Taft-Mojave  net  (Miller 
and  others,  1969b).  Once  again,  the  results  weakly  suggest 
strain  accumulation  rather  than  fault  slippage. 

San  Luis  Obispo-Avenal  Arc 

This  triangulation  arc  beyond  the  northern  edge  of  the 
study  area  (figure  6)  was  established  in  1932  and  reob- 
served  m  1951  and  1962  (Meade,  196S).  Burford  (1966) 
performed  an  extensive  analysis  of  these  data  and  found 


that  they  indicate  right-lateral  fault  slippage  along  the 
San  Andreas  fault. 

Geodimeter  Surveys 

In  1959,  the  California  Department  of  Water  Resources 
started  measuring  a  network  of  15-30  km  length  lines 
along  the  San  Andreas  fault  system  using  a  Model  2A 
Geodimeter  (Hofmann,  1968).  These  data  have  been  dis- 
cussed in  some  detail  by  Hofmann  (1968),  Cherry  and 
Savage  (1972),  Savage  and  Burford  (1973),  Savage  and 
others  (1973),  Greensfelder  and  Bennett  (1973),  and  in 
Appendix  A  of  this  report.  Figure  7  shows  most  of  the 
lines  in  the  Geodimeter  network. 

The  most  pertinent  information  from  these  works  is 
that  strain  is  not  accumulating  along  the  San  Andreas 
fault  north  of  the  vicinity  of  Parkfield  and  that  strain  is 
accumulating  along  the  fault  south  of  Parkfield  (Hof- 
mann, 1968;  Savage  and  Burford,  1973).  The  relative  mo- 
tion between  the  two  sides  of  the  San  Andreas  fault  is 
accommodated  by  fault  creep  at  a  rate  of  between  3  and 
4  cm/year  north  of  Parkfield.  The  Geodimeter  work  indi- 
cates that  fault  creep  ends  about  5  km  south  of  Parkfield. 
Appendix  A  presents  a  dislocation  model  that  explains  the 
Geodimeter  results  in  the  Parkfield  region. 

Fault-Crossing  Survey  Figures 

Beginning  in  1964,  thirty  small  survey  figures  were  es- 
tablished across  faults  of  the  San  Andreas  system,  many 
being  located  where  facilities  of  the  California  Water 
Project  cross  these  faults.  Very  little  analysis  of  the  results 
of  these  figures  has  been  published.  A  compilation  of  the 
results  and  a  brief  analysis  appears  in  Appendix  B.  In  this 
section,  I  will  summarize  the  results  from  Appendix  B. 

The  survey  figures  indicate  right-lateral  movement  on 
the  San  Andreas,  the  Hayward,  and  the  Calaveras  faults 
from  near  Berkeley  to  south  of  Parkfield  with  rates  of 
from  5  to  15  mm/year.  It  is  possible  that  movement  in 
figure  TEM,  south  of  Parkfield,  may  represent  strain  ac- 
cumulation rather  than  fault  creep;  but  there  is  little  doubt 
that  the  movement  in  the  figures  north  of  TEM  represents 
fault  creep. 

From  TEM  south  along  the  San  Andreas  fault  to 
WRIGHT,  the  survey  figures  show  no  movement  larger 
than  the  survey  errors  (a  few  millimeters).  The  figure  at 
DEVIL,  on  the  San  Andreas  fault,  shosvs  no  movement 
until  1970;  but  between  1970  and  1972,  there  is  evidence 
of  right-lateral  movement  at  a  rate  of  several  millimeters 
a  year.  At  the  same  time,  the  two  figures  on  the  San 
Jacinto  fault.  Cf)I.T  and  RIALTO,  indicate  right  lateral 
movement  at  a  rate  of  several  inillimeters  a  year  until  1970 
and  no  movement  after  1970. 

Finally,  two  figures  on  the  Ciarlock  fault  (RANCH- 
TFJON)  indicate  lefi-laleral  movement  at  a  rate  of  8 
mm/year  (Appendix  B).  Sites  on  other  faults  show  no 
movement. 
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Creepmeter  Data 

Since  the  discovery  of  fault  slippage  near  HoUister  in 
1956  (Steinbrugge  and  Zacher,  1960),  much  effort  has 
been  devoted  to  the  study  of  this  phenomenon  (Nason, 
1971,  1973;  Yamashita  and  Burford,  1973).  It  is  now 
known  that  fault  slippage  occurs  on  the  Concord,  the 
Calaveras,  the  Hay  ward,  and  the  San  Andreas  faults  (Na- 
son, 1973).  However,  to  my  knowledge,  there  is  no  pub- 
hshed  evidence  of  slippage  on  any  other  faults. 

The  results  of  the  creepmeter  work  that  are  pertinent  to 
this  study  are  contained  in  Yamashita  and  Burford  (1973) 
and  in  Appendix  A.  At  site  XSCl,  steady  creep  is  occur- 
ring uniformly  at  a  rate  of  about  24  mm/year;  and  at 
XDR 1 ,  steady  creep  occurs  in  steps  at  an  average  rate  of 
1 1  mm/year.  At  XGHl,  the  rate  is  about  2  mm/year;  and, 
at  the  southernmost  site,  XWTl,  there  were  two  creep 
events  between  June  1969  and  February  1973  that  would 
give  an  average  rate  of  about  2  mm/year. 

Thus,  the  data  suggest  a  rather  dramatic  drop  in  the 
creep  rate  between  XDRl  and  XGHl,  a  distance  of  about 
5  km.  However,  it  is  possible  that  the  sites  at  XGHl  and 
XWTl  do  not  cross  ihe  presently  creeping  trace  of  the 
fault.  Two  arguments  suggest  this  may  not  be  a  valid 
explanation.  First,  continual  examination  of  State  High- 
way 46  (by  Division  of  Mines  and  Geology  personnel) 
where  it  crosses  Cholame  Valley  has  not  yielded  evidence 
of  fault  creep  since  the  1966  earthquakes.  If  creep  is  occur- 
ring at  rates  of  1  to  2  cm/year,  it  would  certainly  be 
observed  in  the  highway.  Second,  the  creepmeters  are  in- 
stalled across  the  trace  of  surface  faulting  from  the  1966 
earthquakes,  and  it  seems  reasonable  to  expect  that  creep 
should  occur  on  the  1966  break.  The  Division  of  Mines 
and  Geology  is  conducting  measurements  on  a  survey 
figure  crossing  Cholame  Valley  to  try  to  resolve  this  ques- 
tion. 

Level  Surveys 

The  U.S.  Coast  and  Geodetic  Survey  (now  National 
Ocean  Survey)  has  been  carrying  out  level  surveys  in 
California  since  1906  (Meade  and  Small,  1966).  However, 
very  little  has  been  published  about  the  results  of  the 
leveling.  Parkin  (1948)  discusses  level  changes  around 
Los  Angeles  between  1906  and  1946,  and  he  identifies 
significant  subsidence  near  Long  Beach  and  Santa  Ana. 
Holdahl  (1973)  briefly  considers  the  level  results  at  the 
fault-crossing  survey  sites,  but  he  does  not  present  any  of 
the  data.  Meade  and  Small  (1966)  discuss,  again  briefly, 
some  leveling  results  along  the  San  Andreas  fault. 

Buchanan  and  others  ( 1973)  discuss  elevation  changes 
in  the  western  Transverse  Ranges,  and  they  report  inflec- 
tions in  elevation  changes  that  appear  to  be  associated 
with  faults.  Castle  and  others  (1974)  present  data  from 
repeated  level  surveys  between  western  Los  Angeles  and 
Sandberg  and  Palmdale.  The  results  of  these  repeated  lev- 
els show  that  Sandberg  rose  about  90  mm  relative  to  Los 
Angeles  between  1964  and  1968,  giving  a  rate  of  about  18 


mm/year.  Palmdale  rose  about  180  mm  relative  to  Los 
Angeles  between  1964  and  1968,  giving  a  rate  of  about  36 
mm/year.  They  also  point  out  elevation  differences  that 
they  suggest  might  have  been  precursors  to  the  1971  San 
Fernando  earthquake. 

MATHEMATICAL  MODEL  OF 
THE  SAN  ANDREAS  FAULT 

I  want  to  construct  a  mathematical  model  of  the  San 
Andreas  fault  in  the  Transverse  Ranges  that  incorporates 
the  data  presented  in  "Structural  Characteristics",  "Seis- 
micity",  and  "Geodetic  Characteristics".  The  model  will 
then  be  used  to  attempt  to  account  for  features  of  the 
Quaternary  and  present  tectonic  activity  in  the  Transverse 
Ranges. 

Previous  Work 

Two  earlier  studies  have  approached  the  problem  of 
stress  or  strain  accumulation  in  southern  California. 
Smith  and  Van  de  Lindt  (1969)  calculated  the  changes  in 
strain  energy  due  to  the  introduction  of  stress-free  sur- 
faces (faults)  in  a  pre-stressed  plate  in  plane  strain.  They 
calculated  the  effects  of  movement  on  thirty-one  faults 
from  1812  to  1966,  and  they  found  that  the  strain  energy 
had  increased  in  the  eastern  Transverse  Ranges  (figure  8) . 
They  also  found  an  increase  in  strain  energy  north  of  the 
San  Andreas  fault  near  the  White  Wolf  fault  that  was  only 
partially  relieved  by  the  1952  Arvin-Tehachapi  earth- 
quake. 

This  approach  suffers  from  deficiencies,  several  of 
which  are  discussed  by  the  authors.  The  model  is  two 
dimensional;  hence,  changes  in  strain  cannot  be  calculated 
at  depth,  and  the  orientations  of  potential  failure  planes 
cannot  be  calculated.  Their  model  considers  only  the 
strain  changes  due  to  movements  on  faults;  it  does  not 
consider  the  strain  or  stress  accumulation  occurring 
between  1812  and  1966. 

Rodgers  and  Chinnery  (1973)  attempted  to  deal  with 
the  problem  of  stress  accumulation  since  the  1857  Fort 
Tejon  earthquake.  They  used  a  dislocation  model  of  the 
San  Andreas  fault  that  involved  a  locked  fault  segment 
extending  from  south  of  Parkfield  to  near  San  Bernardino. 
They  found  that  the  maximum  shear  stress  increased  in 
the  Transverse  Ranges  and  in  the  vicinity  of  the  White 
Wolf  fault.  However,  their  work  suffers  from  a  major 
defect.  The  fault  model  they  used  had  a  total  depth  of  100 
kilometers.  A  fault  this  shallow  violates  the  assumption  of 
two  semi-infinite  blocks  moving  past  one  another.  The 
area  where  the  largest  errors  are  introduced  is  exactly  the 
area  of  greatest  interest,  the  area  near  the  bend  in  the  San 
Andreas  fault. 

Review  of  Pertinent  Features 

From  the  discussion  in  the  previous  chapters  and  in  the 
Appendices,  several  important  features  of  the  San  An- 
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Figure  8.  Strain  energy  changes  due 
to  earthquakes,  1812-1967  (from  Smith 
and  Van  de  Lindt,  1969,  p.  1583). 


dreas  fault  emerge  that  must  be  included  in  any  model. 
The  most  obvious  feature  is  the  fact  that  the  fault  has  a 
bend  through  the  Transverse  Ranges.  The  lack  of  creep  on 
the  fault  south  of  the  Parkfield-Cholame  area,  the  lack  of 
seismicity  on  the  San  Andreas  fault  between  Parkfield  and 
San  Bernardino,  and  the  lack  of  movement  seen  in  the 
various  geodetic  networks  must  all  be  included  in  the 
model.  The  high  seismicity  of  the  San  Jacinto  fault  and  the 
relatively  lower  seismicity  of  the  San  Andreas  fault  south 
of  San  Bernardino  must  also  be  considered. 

If  the  mathematical  model  is  a  good  approximation  of 
the  real  world,  the  stress  fields  calculated  from  the  model 
should  account  for  most  of  the  features  of  the  present 
activity  in  the  Transverse  Ranges.  In  particular,  the  east- 
west  trending  thrust  faults  and  the  mechanisms  of  the 
1952  Arvin-Tehachapi  and  the  1971  San  Fernando  earth- 
quakes should  be  reasonably  well  predicted. 


Possible  Mathematical  Models 

The  situation  to  be  modeled  is  that  where  two  relatively 
thin  plates  (the  lithospherc  of  Isacks  and  others,  1968)  are 
moving  past  each  other  along  a  vertical  boundary  (trans- 
form fault),  and  both  plates  are  riding  over  a  much  thicker 
and  weaker  layer  (the  asthenosphere  of  Isacks  and  others, 
1968).  The  model  should  be  three  dimensional  both  to 


accommodate  a  curved  boundary  between  the  plates  and 
to  allow  stresses,  strains,  and  displacements  to  be  calculat- 
ed below  the  surfaces  of  the  plates. 

The  plate  could  contain  a  crack  that  would  approximate 
the  boundary  of  interest.  The  only  closed  solutions  to  this 
problem  of  which  I  am  aware  involve  single  straight 
cracks  where  the  plate  is  assumed  to  be  in  a  state  of  plane 
strain.  Thus,  a  curved  boundary  could  not  be  modeled; 
there  would  be  no  variation  of  the  calculated  fields  with 
depth,  and  inhomogeneities  could  not  be  included. 

Another  class  of  models  is  that  involving  cracks  or 
dislocations  in  half-spaces.  The  equations  for  displace- 
ments, strains,  and  stresses  due  to  a  crack  or  dislocation 
in  a  homogeneous  half-space  have  been  developed  in 
closed  form  and  thus  could  be  used  fairly  easily.  Also, 
since  dislocation  models  in  particular  are  reasonably  good 
representations  of  faults  (Chinnery,  1961,  1963),  models 
involving  dislocations  are  desirable.  The  major  drawback 
to  a  dislocation  model  is  that  the  equations  for  a  layered 
medium  exist  only  in  integral  form,  and  the  solutions 
require  numerical  integrations  with  their  attendant  prob- 
lems of  cost  and  accuracy.  On  the  other  hand,  a  disloca- 
tion model  would  allow  a  curved  boundary  between  two 
blocks  to  be  modeled;  and  it  would  allow  the  displace- 
ments, strains,  and  stresses  to  be  calculated  below  the 
surface. 
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Some  work  has  been  done  on  the  two-dimensional 
problem  of  a  dislocation  in  a  layered  medium  (McGiniey, 
1969;  Chinnery  and  Jovanovich,  1972;  Jovanovich  and 
others,  1974a,b).  This  work  could  be  used  to  estimate  the 
effects  of  layering  on  the  results  from  a  three-dimensional 
homogeneous  half-space  model. 

This  is,  then,  the  approach  I  will  take.  A  model  of  the 
San  Andreas  fault  in  southern  California  will  be  con- 
structed using  a  dislocation  surface  in  a  homogeneous 
half-space,  and  the  effect  of  layering  on  the  results  will  be 
estimated  from  the  two-dimensional  models. 


Details  of  the 
Selected  Dislocation  Model 

Dislocation  theory  has  been  used  for  straight  finite 
fauhs  (Chinnery,  1961,  1963),  for  finite  faults  with  bends 
(Rodgers  and  Chinnery,  1973),  and  for  infinitely  long 
straight  faults  (Chinnery,  1961).  In  the  case  of  a  fault 
with  a  bend  in  it,  the  fault  should  be  very  long  to  remove 
the  effects  of  the  ends  from  the  area  of  interest.  I  would 
like  to  use  the  expressions  developed  by  Chinnery  (1961, 
1963)  to  model  the  present  situation,  because  these  ex- 
pressions are  closed  and  fairly  simple.  Unfortunately,  to 
use  these  expressions,  certain  assumptions  are  immediate- 
ly required.  Chinnery's  expressions  were  developed  for  a 
straight,  finite,  vertical  strike-slip  fault  in  an  isotropic, 
homogeneous,  linearly  elastic  half-space. 

The  expressions  for  a  straight  strike-slip  fault  can  be 
used  for  the  case  of  a  curved  fault  if  we  approximate  the 
fault  with  a  series  of  straight  fault  segments.  The  displace- 
ments, strains,  and  stresses  at  an  observation  point  can 
then  be  determined  by  summing  the  fields  at  the  observa- 
tion point  that  are  produced  by  each  of  the  straight  seg- 
ments. Care  must  be  taken,  of  course,  to  ensure  that  all  the 
fields  are  transformed  into  the  same  coordinate  system. 
There  will  be  some  perturbations  introduced  into  the  cal- 
culated fields  due  to  the  use  of  finite  segments;  but,  if 
enough  segments  are  used  and  the  grid  of  observation 
points  is  spaced  properly,  this  effect  can  be  greatly  re- 
duced. 

Details  of  the  Half-Space 

The  Transverse  Ranges  province  is  anything  but  an 
isotropic,  homogeneous,  linearly  elastic  half-space.  The 
presence  of  mountains,  faults,  very  different  rock  types, 
and  deep  subhorizontal  discontinuities  must  change  the 
calculated  elastic  fields.  McGiniey  (1969)  and  Ben- 
Menahem  and  others  (1970)  have  shown  that,  at  dis- 
tances of  less  than  about  20°,  the  half-space  adequately 
represents  the  curved  earth.  Since  the  area  I  am  interested 
in  is  less  than  3°  across,  at  least  the  assumption  of  a  flat 
earth  should  not  be  a  problem. 

The  problem  of  inhomogeneities  is  not  as  easy  to  dis- 
miss. Lateral  changes  in  rock  types  will  affect  the  calculat- 
ed fields,  but  presumably  the  results  would  be  fairly  local 


in  extent.  These  lateral  inhomogeneities  will  almost  cer- 
tainly lead  to  local  stress  concentrations  near  the  surface; 
but  if  their  vertical  extents  are  not  great,  their  effects 
should  be  confined  to  areas  near  the  inhomogeneities.  In 
fact,  the  local  stress  concentrations  from  these  in- 
homogeneities probably  determine  the  local  areas  most 
susceptible  to  earthquakes. 

Potentially  more  important  is  the  problem  presented  by 
vertical  inhomogeneities.  The  plate  tectonics  hypothesis  as 
advanced  by  Isacks  and  others  (1968)  divides  the  outer 
part  of  the  earth  into  three  zones  of  very  different  physical 
properties.  The  upper  zone,  the  lithosphere,  is  approxi- 
mately 100  km  thick,  is  fairly  "strong",  and  exhibits  more 
or  less  elastic  behavior.  The  middle  zone,  the  astheno- 
sphere,  is  relatively  much  thicker  and  "weaker"  and 
behaves  more  or  less  plastically.  The  third  zone,  the 
mesosphere,  is  relatively  "stronger"  than  the  astheno- 
sphere  and  may  or  may  not  exhibit  viscous  properties. 

On  a  finer  scale,  the  upper  30  or  so  kilometers  of  the 
lithosphere  in  southern  California  is  also  divided  into  two 
or  more  layers  or  zones  (Press,  1960;  Healy,  1963).  For 
example,  the  crustal  model  derived  by  Press  (1960)  has 
two  layers  in  the  crust  above  the  crust-mantle  boundary 
(the  Moho).  East  of  the  San  Andreas  fault,  the  Moho  is 
about  50  km  deep  (Press,  1960);  while  west  of  the  fault, 
the  Moho  is  only  about  30  km  deep  (Healy,  1963). 

Jovanovich  and  others  (1974b)  have  considered  the 
effects  of  layers  on  the  far  field  strains  and  displacements 
from  a  dislocation.  They  used  a  two-dimensional  model 
of  a  point  source  in  a  multi-layered  medium  and  numeri- 
cally evaluated  the  strain  and  displacement  fields.  Their 
results  showed  that  a  thick  "weak"  layer  underlying  a 
thick  "strong"  layer  increases  by  a  factor  of  2  or  less  the 
strain  and  tilt  fields  one  would  get  by  using  an  identical 
source  in  a  homogeneous  half-space.  They  also  showed 
that  a  relatively  thin  "weak"  layer  between  two  "strong- 
er" layers  drastically  alters  both  the  sign  and  magnitude 
of  the  strains  calculated  from  an  identical  source  in  a 
homogeneous  half-space. 

Thus,  based  on  the  previous  work,  the  effect  of  the 
lithosphere  on  the  elastostatic  fields  calculated  from  my 
model  would  be  to  increase  the  calculated  fields  by  per- 
haps as  much  as  a  factor  of  2.  The  effect  of  a  thinner, 
"weaker"  layer  on  the  calculated  fields  is  somewhat  un- 
predictable but  could  be  drastic.  There  is  some  weak  evi- 
dence for  low- velocity  zones  in  the  crust  (Landisman  and 
others,  1971),  but  the  presence  of  such  zones  is  still  far 
from  certain  (Healy,  1971;  Warren  and  Healy,  1973). 
Therefore,  all  that  can  be  said  is  that  a  thin  low-velocity 
("weak")  layer  might  exist  under  the  study  area;  and  if 
it  exists,  its  effect  on  the  calculated  fields  is  unpredictable. 
However,  I  do  not  feel  that  the  present  evidence  for  a 
crustal  low-velocity  zone  is  very  good.  For  these  reasons, 
the  large-scale  features  of  the  elastostatic  fields  calculated 
from  the  homogeneous  model  are  probably  reasonable 
approximations  of  the  fields  that  exist  in  the  earth. 
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The  basic  model,  then,  is  one  of  a  curved  dislocation 
surface  in  an  isotropic,  homogeneous,  linearly  elastic  half- 
space.  Recall,  however,  that  I  am  trying  to  model  the  case 
of  two  plates  moving  past  one  another  along  a  curved 
boundary.  For  the  half-space  model  to  fit  this  situation, 
it  seems  reasonable  to  require  that  the  boundary  extend 
vertically  to  infinity.  Thus,  we  now  have  two  quarter- 
spaces  moving  past  one  another  on  a  curved  boundary. 

Details  of  the  Boundary 

Let  us  now  consider  the  details  of  the  boundary  between 
the  two  blocks.  First  of  all,  I  am  really  only  interested  in 
the  small  region  in  the  vicinity  of  the  curve  in  the  bound- 
ary, so  it  seems  reasonable  to  place  the  two  ends  of  the 
boundary  at  a  great  distance  from  the  bend.  This  will 
completely  remove  the  end  effects  of  the  dislocation  from 
the  effects  due  to  the  bend. 
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Figure  9.        Dislocation  surface  and  its  fit  to  faults. 


The  observed  geological,  seismological,  and  geodetic 
data  suggest  that,  at  the  present  time,  the  section  of  the 
San  Andreas  fault  containing  the  bend  is  not  moving  but 
that  the  straight  sections  on  either  side  of  the  bend  are 
moving.  These  data  led  Allen  (1968)  to  state  that  this  part 
of  the  fault  is  "locked".  What,  then,  is  the  mechanism  of 
"locking"?  Intuitively,  one  would  expect  the  blocks  to 
experience  greater  difficulty  in  moving  around  a  bend 
than  in  moving  along  a  straight  boundary.  Thus,  when 
strain  accumulation  at  either  end  of  the  bend  has  reached 
a  high  enough  level,  the  "locked"  section  of  the  fault  will 
fail  and  produce  quite  a  large  earthquake.  Presumably,  the 
amount  of  movement  in  the  earthquake  will  relieve  most 
or  all  of  the  accumulated  strain. 

The  slip  on  the  moving  sections  of  the  fault  on  either 
side  of  the  bend  is  probably  occurring  as  very  slow  visco- 
elastic  creep  with  no  release  of  seismic  energy  below  about 
20  km  and  as  intermingled  creep  and  small  earthquakes 
above  20  km  (Scholz  and  others,  1969;  Brune,  1974).  The 
slip  behavior  above  20  km  is  similar  to  the  stick-slip 
behavior  of  rock  specimens  in  the  laboratory  (Byerlee  and 
Brace,  1 968 ) .  Thus,  it  seems  reasonable  to  assume  that  the 
"locking"  mechanism  around  the  bend  is  probably  a  large 
-scale  example  of  stick-slip  behavior. 

The  depth  of  the  "locked"  zone  is  also  important.  It 
seems  reasonable  to  assume  that  the  depth  of  the  stick-slip 
behavior  around  the  bend  is  about  the  same  (20  km)  as 
that  of  the  straight  sections.  Additionally,  it  seems  reason- 
able to  say  that  the  viscoelastic  creep  under  the  straight 
fault  sections  continues  under  the  "locked"  section. 

The  complete  model  is  shown  in  figure  9,  and  figure  10 
shows  the  coordinate  system  used.  Note  that  most  of  the 
fault  surface  is  a  zone  of  stable  sliding  and  that  the  locked 
zone  is  shallow  compared  to  the  total  fault  depth.  In  the 
next  section,  the  parameters  used  for  the  model  will  be 
discussed  together  with  the  fit  of  the  model  to  the  real 
world. 
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Parameters  of  the 
Dislocation  Model 

To  model  the  situation  described  previously,  I  need  to 
determine  as  precisely  as  possible  the  extent  of  the  active 
fault  surface  between  the  two  blocks,  the  ends  of  the 
locked  segment,  the  depth  of  the  locked  segment,  and  the 
amount  of  displacement  on  the  fault.  The  lengths  of  the 
end  segments  and  the  depth  of  the  bottom  of  the  moving 
zone  will  be  made  very  large  as  discussed  in  the  last  sec- 
tion. This  will  remove  the  effects  of  the  ends  and  bottom 
of  the  fault.  The  other  parameters  can  be  determined  from 
the  geological,  seismological,  and  geodetic  data  discussed 
earlier. 

The  boundary  between  the  two  blocks  is  distinct  from 
the  vicinity  of  San  Juan  Baulista  to  near  Cajon  Pass.  All 
of  the  evidence  indicates  that  the  San  Andreas  fault  is  the 
only  through-going,  continuously  or  nearly  continuously 
active  fault;  and  it  is  reasonable  to  say  that  it  represents 
the  boundary  between  the  blocks.  However,  south  of  Ca- 
jon Pass,  the  boundary  is  not  so  obvious.  Allen  ( 1957)  was 
unable  to  trace  the  San  Andreas  fault  through  San  Gor- 
gonio  Pass. 
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The  San  Andreas  fault  is  not  very  active  seismically 
between  Parkfield  and  San  Gorgonio  Pass,  and  there  is 
very  little  evidence  of  movement  along  the  fault.  However, 
the  San  Jacinto  fault  has  had  at  least  five  earthquakes  with 
magnitudes  greater  than  6.0  associated  with  it  (figure  3), 
there  has  been  surface  breakage  along  most  of  its  length 
(figure  4),  and  there  is  evidence  of  continuing  movement 
along  the  fault  (Appendix  B).  For  these  reasons,  it  seems 
reasonable  to  assume  that  the  San  Jacinto  fault  south  of 
Cajon  Pass  is  the  boundary  between  the  blocks  and  that 
there  is  some  complex  interaction  between  the  San  An- 
dreas and  the  San  Jacinto  faults  near  Cajon  Pass. 

I  have  fitted  straight  fault  segments  to  the  traces  of  the 
San  Andreas  and  San  Jacinto  faults  as  shown  on  the 
1:250,(XX)  sheets  of  the  Division  of  Mines  and  Geology 
Geologic  Map  of  California.  Figure  9  shows  the  traces  of 
the  faults  and  the  model  fit  to  the  faults.  The  agreement 
between  the  model  and  the  mapped  faults  is  very  good 
from  Parkfield  to  south  of  San  Bernardino. 

The  northern  end  of  the  locked  segment  is  probably 
determined  to  within  10  km  by  the  seismicity,  the  abrupt 
cessation  of  fault  creep,  and  geodetic  evidence  (Appendix 
A).  The  northern  end  of  the  locked  segment  was  placed 
near  Parkfield.  The  southern  end  is  less  well  determined. 
There  is  no  published  evidence  of  surface  creep  along  the 
San  Jacinto  fault,  but  there  does  appear  to  have  been 
surface  rupture  associated  with  the  large  earthquakes 
along  the  fault  trace  south  of  San  Bernardino.  The  south 
end  of  the  locked  zone  was  placed  near  San  Bernardino, 
and  it  is  probably  accurate  to  within  30  km.  It  should  be 
noted  that  the  exact  location  of  this  end  of  the  locked 
segment  is  less  critical  than  the  location  of  the  north  end, 
since  I  am  mainly  interested  in  the  stress  accumulation  in 
the  western  Transverse  Ranges  province.  The  effects  of 
changing  the  location  of  the  southern  end  of  the  segment 
will  not  be  large  in  the  area  of  interest. 

The  depth  of  the  bottom  of  the  locked  zone  is  more 
difficult  to  determine.  It  seems  reasonable  to  assume  that 
the  depth  of  the  locked  segment  corresponds  roughly  to 
the  depth  of  the  fault  segment  that  exhibits  stick-slip 
behavior  in  central  California.  The  vertical  extent  of  this 
segment  is  well  defined  by  microearthquakes  (Brune  and 
Allen,  1967;  Bolt  and  Miller,  1971 ),  and  the  bottom  of  the 
zone  seems  to  be  between  12  and  15  km.  The  aftershocks 
of  the  1966  Parkfield  earthquake  were  no  deeper  than  15 
km  (Eaton  and  others,  1970),  and  a  few  microearth- 
quakes located  near  the  HUGHES  fault-crossing  survey 
site  (figure  B-1)  were  no  deeper  than  15  km  (Brune  and 
Allen,  1967).  Thus,  a  depth  of  15  km  for  the  bottom  of 
the  locked  segment  of  the  fault  seems  to  be  indicated  by 
the  seismicity  data.  The  bottom  of  the  locked  segment 
becomes  progressively  shallower  to  either  end  to  smooth 
out  end  effects  to  some  degree.  Appendix  A  contains  a 
discussion  of  the  degree  of  control  on  the  depth  of  the 
bottom  of  the  locked  segment  near  Parkfield.  Figure  9 
shows  the  shape  of  the  locked  segment  used  in  this  study. 


Due  to  the  uncertainty  of  the  depth  of  the  locked  segment, 
results  from  models  with  depths  of  15  and  30  km  will  be 
used. 

The  amount  of  displacement  on  the  moving  segments  is 
fairly  well  determined  from  the  geodetic  work.  The  es- 
timated displacement  rate  on  the  San  Andreas  fault  north 
of  Parkfield  is  about  3  cm/year  (Savage  and  Burford, 
1973).  The  estimated  displacement  rate  across  the  Impe- 
rial Valley  may  be  as  much  as  8  cm/year  (Whitten,  1960; 
Scholzand  Fitch,  1969,  1971;  Savage  and  Burford,  1971). 
Atwater  and  Molnar  (1973)  determined  a  post^.5-mil- 
lion-year  rate  of  slip  of  5.5  cm/year  from  sea-fioor 
spreading  data.  For  the  present  study,  a  rate  of  5  cm/year 
is  used. 

It  should  be  instructive  to  look  at  the  stress  accumula- 
tion after  the  1857  earthquake  and  before  the  1952  Arvin- 
Tehachapi  earthquakes;  thus,  the  displacement  on  the 
moving  segments  for  the  94  years  1858-1952  at  a  rate  of 
5  cm/yr  is  4.7  meters.  This  is  the  value  of  the  displacement 
used  in  this  study.  It  should  be  noted  that  the  value  of  the 
displacement  does  not  affect  the  shape  of  the  calculated 
displacement,  strain,  and  stress  fields.  It  affects  only  the 
value  of  these  fields.  Table  1  gives  the  parameters  of  the 
two  dislocation  models. 

For  the  purposes  of  this  calculation,  Lame's  constants 
are  set  equal  (i.e.,  Poisson's  Ratio  =  0.25)  and  have  a 
value  of  3.10"  dynes/cm'.  This  is  a  fairly  common  as- 
sumption that  simplifies  the  calculations  (Chinnery,  1961, 
1963;  Mansinha  and  Smylie,  1971). 

The  displacement  on  each  fault  segment  was  deter- 
mined by  Ui  =  U/cos  0i,  where  U  =  4.7  meters,  u,  =  the 
displacement  on  the  i""  segment,  and  0,  is  the  angle  the  i'*" 
segment  makes  with  the  Y,  axis.  This  form  of  the  displace- 
ment is  required  by  the  assumption  of  constant  displace- 
ment in  the  blocks  at  some  distance  from  the  fault.  If,  say, 
the  displacement  were  the  same  on  each  fault  segment, 
this  would  mean  the  displacement  parallel  to  the  Yj  axis 
would  have  to  vary. 

The  displacement  on  each  segment  was  constant  over 
the  whole  area  of  the  segment.  While  this  does  produce 
discontinuities  in  the  calculated  fields  at  the  edges  of  each 
segment,  it  is  felt  that  this  is  not  a  problem  so  long  as  the 
observation  points  are  greater  than  10  km  from  the  fault 
trace.  This  introduces  some  short-wavelength  changes  in 
the  calculated  fields,  but  proper  grid  spacing  can  eliminate 
these. 

CALCULATED  VERTICAL 
DISPLACEMENT  FIELDS 

Figures  11  and  12  show  the  vertical  displacements  of 
the  ground  surface  produced  by  models  1  and  2.  The 
vertical  displacement  fields  from  each  model  have  the 
same  overall  shape,  although  there  are  some  minor  local 
variations.  There  is  general  uplifi  over  much  of  the  study 
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Toble    1.      Parameters  of  the  models. 


Segment 

Sefri'ifnt 

length 

Displacement 

nui7ib<^r 

2L  (km) 

V  (m) 

0 

100,000 

4  70 

1 

50 

4.70 

2 

SO 

4.70 

3 

3.62 

4.75 

4 

9.50 

4.70 

5 

4.74 

4.71 

6 

612 

4.72 

7 

5.74 

4.70 

8 

6,38 

4.70 

9 

3.12 

4.72 

10 

8.00 

4.70 

U 

1.05 

4.76 

12 

8.38 

4.71 

13 

7.38 

4.75 

14 

8.00 

4.80 

15 

3.88 

4.75 

16 

4.24 

4.89 

17 

4.24 

5.00 

18 

13.38 

557 

19 

9.74 

6.43 

20 

11.74 

5.60 

21 

4.24 

7.24 

22 

3.50 

6.01 

23 

9.00 

5.32 

24 

32.74 

5.51 

25 

15.24 

5.30 

26 

9.00 

532 

27 

7.74 

5.27 

28 

42.24 

5.14 

29 

42.86 

4.82 

30 

100,000 

4.76 

Depth  to 

Depth  to  top 

bottom 

dl  (km) 

d2  (km) 

D  (km) 

(1 

0 

1(X).(X)0 

5 

10 

lOO.IXX) 

10 

20 

1(XI.(XI0 

15 

30 

100.000 

IS 

30 

100.000 

15 

30 

100.000 

15 

30 

100,000 

15 

30 

100.000 

15 

30 

100,000 

15 

30 

100.000 

15 

30 

1IXI.(X)0 

15 

30 

luo.ooo 

15 

30 

100,000 

IS 

30 

1(X),000 

15 

30 

100,000 

15 

30 

lOO.lXX) 

15 

30 

100,000 

15 

30 

100,000 

15 

30 

100,000 

15 

30 

KXI.IXX) 

15 

30 

i(X),(xx) 

15 

30 

UXI,(XK) 

15 

30 

100.000 

15 

30 

100,000 

15 

30 

i(X).(xx) 

15 

30 

1(K),(XK) 

15 

30 

1(X),(XK) 

15 

30 

1IX),(XX) 

10 

20 

i(x).ooo 

5 

10 

UX).(KX) 

0 

0 

i(xi.(xx) 

area,  with  the  largest  uplift  occurring  in  the  vicinity  of 
Mount  Pinos  and  Frazier  Mountain  with  a  magnitude  of 
between  40  and  45  cm.  This  suggests  a  maximum  rate  of 
uplift  from  1858  to  1952  of  about  4.3  to  4.8  mm/yr  in  the 
vicinity  of  Mount  Pinos.  The  relative  uplift  is  generally 
much  lower  near  the  Carrizo  Plain,  in  the  vicinity  of  San 
Bernardino,  and  in  the  vicinity  of  Lancaster  and  Palm- 
dale. 

Two  points  should  be  made  about  the  vertical  displace- 
ments. First,  the  area  of  greatest  relative  uplift  is  crudely 
reflected  in  the  topography.  Second,  the  displacement 
fields  from  each  model  resemble  one  another  much  more 
than  one  would  expect. 

The  obvious  question  is  whether  the  rather  large  values 
of  the  calculated  vertical  displacements  are  due  to  the  use 
of  a  homogeneous  half-space.  Jovanovich  and  others 
(1974)  have  shown  that,  for  a  buried  point  source,  the 
vertical  displacements  close  to  the  source  are  smaller  for 
the  half-space  model  than  for  a  model  with  a  thin  weak 
layer  at  the  surface  and  that  this  half-space  model  has 
about  the  same  value  as  a  model  with  a  buried  thin  weak 
layer.  At  greater  distances  from  the  source,  the  model  with 
a  buried  thin  weak  layer  produces  larger  vertical  displace- 
ments than  the  half-space  model  or  the  model  with  a  thin 
surface  layer.  Thus,  in  both  layered  cases,  the  homogene- 
ous half-space  produces  a  minimum  estimate  of  the  verti- 
cal displacement. 


If  the  calculated  vertical  displacements  are  a  reasonable 
representation  of  the  "real"  earth,  it  is  necessary  to  ask 
whether  some  mechanism  can  operate  to  remove  them  in 
100  years  or  less.  Three  such  mechanisms  come  to  mind: 
viscoelastic  relaxation,  isostatic  adjustment,  and  erosion. 

Singh  and  Rosenman  (1974)  have  shown  that,  using 
viscoelastic  models  of  the  earth,  displacements  do  not  tend 
to  die  out  with  time.  Therefore,  I  would  not  expect  that 
the  vertical  deformation  produced  by  either  model  1  or 
model  2  would  decrease  in  time  due  to  a  viscoelastic  relax- 
ation process.  The  only  other  ways  to  remove  the  vertical 
changes  would  be  through  either  erosion  or  isostatic  ad- 
justment. 

Hanks  (1974)  states  that  there  is  no  evidence  for  a  root 
under  the  San  Gabriel  Mountains  (figures  11  and  12),  and 
this  suggests  that  isostatic  adjustment  is  not  occurring. 
His  discussion  does  not  rule  out  the  possibility  of  some 
sort  of  broad  regional  compensation  such  as  "an  elastic 
layer  deformed  by  the  vertical  load  exerted  by  the  uplifted 
mass  overlying  a  fluid  substratum."  (Hanks,  1974,  p. 
1226).  However,  the  stresses  produced  from  this  model 
are  not  consistent  with  the  observed  north-south  com- 
pression in  the  Transverse  Ranges  (Hanks,  1974),  and 
Hanks  suggests  that  the  uplifted  mass  may  be  supported 
by  the  horizontal  compression.  He  then  proceeds  to  derive 
a  value  of  about  1  kb  for  the  mean  horizontal  stress  neces- 
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sary  to  hold  up  the  elevated  mass.  But  this  estimate  of  the 
horizontal  stress  is  about  100  times  the  horizontal  stresses 
calculated  from  my  model.  The  half-space  model  may 
underestimate  the  horizontal  stresses,  but  it  does  not  seem 
reasonable  to  say  that  the  stresses  are  underestimated  by 
2  orders  of  magnitude.  Thus,  there  may  be  some  regional 
isostatic  compensation  operating  to  help  to  account  for 
the  difference  between  Hanks'  (1974)  estimate  of  the 
horizontal  stress  and  the  vertical  horizontal  stresses  cal- 
culated from  my  model. 

The  elevation  increases  could  also  be  removed  by  ero- 
sion if  denudation  is  taking  place  at  a  rate  of  about  4  to 
5  mm/yr.  However,  the  highest  observed  denudation  rate 
in  the  United  States  is  only  0.8  mm/yr  in  the  Eel  River 
Basin  in  the  northern  Coast  Ranges  of  California  (Judson 
and  Ritter,  1964).  Since  neither  isostatic  adjustment  nor 
erosion  seems  to  be  occurring  at  a  rate  fast  enough  to 
remove  the  elevation  increases  indicated  by  the  model,  we 
are  left  with  the  possibility  that  the  general  form  of  the 
uplift  of  the  Tehachapi  Mountains  and  Transverse  Ranges 
might  be  produced  by  the  bend  in  the  San  Andreas  fault. 

CALCULATED  STRESS  FIELDS 

All  six  stress  components  at  each  observation  point 
were  calculated  for  the  two  models.  Due  to  the  difficulty 
of  attaching  physical  meaning  to  individual  components 
of  the  stress  tensor,  the  individual  components  will  not  be 
presented.  Instead,  other  representations  that,  in  a  sense, 
summarize  the  stress  tensor  will  be  used.  Isochromatics, 
or  curves  of  constant  maximum  shear  stress,  will  be  pre- 
sented. The  isochromatics  can  be  used  to  determine  areas 
of  potential  failure.  Isopachs,  or  curves  of  constant  mean 
stress,  will  be  presented  because  they  can  be  used  to  help 
determine  areas  of  potential  failure  (Chinnery,  1966a,b). 
To  determine  directions  of  potential  failure,  stress-trajec- 
tories (isostatics),  or  curves  whose  directions  at  any  point 
are  the  directions  of  the  principal  stress  axes,  will  be  pre- 
sented at  the  surface.  Also,  slip-lines,  a  set  of  orthogonal 
curves  whose  directions  at  any  point  represent  the  direc- 
tions of  the  planes  across  which  the  maximum  shear 
stresses  act,  will  be  presented  for  the  models  at  the  surface. 
This  nomenclature  follows  that  of  Jaeger  and  Cook  (1969, 
p.  17). 

In  general,  it  is  not  practical  to  present  stress-trajecto- 
ries and  slip-planes  at  depth  because  the  principal  stress 
axes  do  not  lie  in  one  plane  and  the  planes  of  maximum 
shear  stress  do  not  intersect  one  plane  at  90°.  Both  of  these 
conditions  hold  only  at  a  free  surface.  I  will  not  attempt 
to  present  here  principal  stress  axes  or  potential  failure 
planes  below  the  surface. 

Strictly  speaking,  the  calculated  stress  fields  presented 
here  are  changes  in  some  (presumably  non-zero)  pre- 
existing stress  field.  However,  I  feel  that  the  calculated 
stress  fields  may  represent  the  present  stress  field  reason- 
ably well.  Rosenman  and  Singh  (1973)  have  shown  that 
all  the  components  of  the  stress  tensor  will  fall  ofT  expo- 


nentially in  time  if  faulting  occurs  in  a  viscoelastic  half- 
space.  For  example,  if  the  relaxation  time  is  50  years,  in 
100  years  the  stresses  will  have  fallen  to  about  10%  of  the 
values  at  the  time  of  failure;  and  if  the  relaxation  time  is 
100  years,  the  values  will  have  fallen  to  between  40%  and 
50%  of  the  initial  stress.  This  suggests  that  any  stresses 
not  relieved  by  the  1857  earthquake  may  have  fallen  to  a 
relatively  low  value  by  1952.  Therefore,  the  stress  field 
primarily  responsible  for  the  1952  Arvin-Tehachapi 
earthquakes  would  have  been  the  stress  field  due  to  the 
lack  of  movement  on  the  locked  segment  of  the  San  An- 
dreas fault. 

Results  from  Model  1 

The  stress  data  calculated  from  model  1  will  be  present- 
ed first.  The  stresses  were  calculated  for  depths  of  0,10, 
and  20  km.  A  regional  stress  field  of  the  form  cr,  =  cti 
=  CTj  =  cr  =    pgy„  where: 

p  is  the  density  of  the  material, 

g  is  the  acceleration  of  gravity,  and 

o",  are  the  principal  stresses  of  the  regional  field, 

should  be  added  to  the  results. 

This  stress  field  is  the  lithostatic  pressure  at  a  depth  y, 
and  represents  uniform  compression.  By  leaving  this  stress 
field  out,  we  can  study  the  effects  of  the  stress  deviations 
produced  by  the  bend  that  are  related  to  failure  criteria 
(Jaeger  and  Cook,  1969,  p.  33).  Note  that  such  a  regional 
lithostatic  pressure  will  not  affect  the  orientation  of  the 
principal  stress  axes.  Of  the  quantities  being  presented, 
only  the  values  of  the  average  stress  are  affected  by  the 
lithostatic  pressure. 

Figure  13  shows  the  isochromatics  (lines  of  constant 
maximum  shear  stress)  calculated  from  model  1  at  the 
ground  surface.  The  maximum  shear  stress  is  quite  large 
near  the  fault  trace;  and  in  some  places,  the  values  are  very 
much  larger  than  the  13  bars  (13  x  10' dynes/cm')  shown. 
Note  particularly  the  lobe  to  the  north  of  the  White  Wolf 
fault  and  the  large,  diffuse  lobe  including  the  Newport- 
Inglewood  fault  zone. 

Figure  14  shows  contours  of  constant  values  of  the  sum 
of  the  principal  stresses.  This  quantity  is  three  times  the 
change  in  mean  stress  (Jaeger  and  Cook,  1969,  p.  17). 
This  change  in  mean  stress  shows  areas  where  the  litho- 
static compression  discussed  earlier  is  either  increased 
(positive  values)  or  decreased  (negative  values).  There  is 
a  broad  zone  of  increased  compression  running  north- 
south  across  figure  14  with  areas  of  decreased  compres- 
sion near  the  Big  Pine  fault. 

Figure  1 5  shows  the  calculated  stress  trajectories  at  the 
ground  surface.  The  axis  of  greatest  compression  (dashed 
lines)  trends  generally  north-south  throughout  most  of 
the  area,  but  this  axis  acquires  a  more  southeast-to-east 
trend  between  the  Newpwrt-Inglewood  fault  zone  and  the 
San  Jacinto  fault  zone. 
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Figure  16  shows  surface  traces  of  the  slip-planes 
(planes  across  which  the  maximum  shear  stress  acts)  at 
the  surface.  These  planes  make  angles  of  45°  with  the  axis 
of  greatest  compression,  and  they  are  the  planes  on  which 
failure  should  theoretically  occur.  However,  laboratory 
experiments  suggest  that  failure  planes  typically  occur  at 
angles  of  between  30°  and  45°  to  the  principal  compression 
axis  (Jaeger  and  Cook,  1969).  Note  that  the  right-lateral 
slip  lines  in  figure  16  tend  to  be  parallel  or  sub-parallel  to 
most  of  the  right-lateral  faults  in  the  area  such  as  the 
Newport-Inglewood  zone,  the  Whittier  fault,  and  the  San 
Gabriel  fault.  The  left-lateral  slip  lines  are  sub-parallel  to 
the  Garlock,  the  Big  Pine,  the  Pleito,  and  the  Frazier 
faults. 

Figure  17  shows  the  isochromatics  calculated  from 
model  1  at  a  depth  of  10  km.  The  contours  are  similar  to 
those  of  figure  13.  Note  the  lobe  near  the  western  end  of 
the  White  Wolf  fault  and  the  relatively  steep  stress  gradi- 
ent parallel  to  and  on  either  side  of  the  San  Andreas  fault. 

Figure  1 8  shows  contours  of  equal  values  of  the  sum  of 
the  principal  stresses  at  a  depth  of  10  km.  Figure  18  resem- 
bles figure  14  except  that  there  is  no  area  of  decreased 
compression  near  the  Big  Pine  fault. 

Figure  19  shows  the  isochromatics  at  a  depth  of  20  km 
calculated  from  model  1.  Note  the  similarity  to  figures  13 
and  15.  These  figures  indicate  that  the  maximum  shear 
stress  is  fairly  constant  with  depth. 

Figure  20  shows  isopachs  calculated  from  model  1  at  a 
depth  of  20  km.  Once  again,  the  features  are  generally 
similar  to  those  of  figures  14  and  18,  but  there  are  some 
local  variations.  Note,  for  example,  the  change  in  location 
of  the  zero  contour  near  the  Big  Pine-San  Andreas  junc- 
tion. 

Results  from  Model  2 

Model  2  has  the  depth  of  the  bottom  of  the  locked  zone 
set  at  30  km.  The  other  parameters  are  the  same  as  for 
model  1.  Figure  21  shows  the  isochromatics  calculated 
from  model  2  at  the  ground  surface.  Comparing  figure  21 
with  figure  13,  we  can  see  that  the  3-bar  contour  is  some- 
what farther  from  the  San  Andreas  fault  for  model  2,  that 
the  gradients  in  the  stress  field  are  less  for  model  2,  and 
that  the  maximum  shear  stresses  are  lower  near  the  fault 
for  model  2.  The  general  form  of  the  isochromatics  is 
about  the  same,  with  lobes  of  increased  maximum  shear 
stress  developed  north  of  the  White  Wolf  and  south  of  the 
Newport-Inglewood  faults,  and  lows  extending  parallel  to 
the  Malibu-Santa  Monica  zone  and  the  White  Wolf  fault. 

Figure  22  shows  the  isopachs  calculated  from  model  2 
at  the  ground  surface.  There  is  an  area  of  decreased  pres- 
sure running  across  the  figure  parallel  to  the  White  Wolf 
and  the  Garlock  faults.  There  are  also  areas  of  decreased 
pressure  to  the  north  of  the  San  Andreas  fault  in  the 


northeast  part  of  the  figure  and  a  smaller  area  in  the 
northwest  corner  of  the  figure. 

Figure  23  shows  the  stress  trajectories  derived  from 
model  2  at  the  ground  surface.  Note  that,  over  most  of  the 
area  covered  by  figure  23,  the  principal  compressive  stress 
is  oriented  north-south.  The  stress  trajectories  derived 
from  model  2  are  much  smoother  than  those  derived  from 
model  1. 

Figure  24  shows  the  slip  lines  calculated  from  model  2 
at  the  ground  surface.  The  right-lateral  slip  lines  are  gen- 
erally sub-parallel  to  the  San  Gabriel,  the  Whittier,  the 
Chino,  and  the  Newport-Inglewood  faults,  while  the  left- 
lateral  slip  lines  tend  to  make  angles  of  20  to  30°  with  the 
left-lateral  faults,  such  as  the  White  Wolf,  the  Garlock, 
the  Santa  Ynez,  and  the  Malibu-Santa  Monica  faults.  Part 
of  the  Big  Pine  fault  is  sub-parellel  to  a  left-lateral  slip 
line.  The  slip  lines  derived  from  model  2  are  somewhat 
smoother  than  those  derived  from  model  1. 

Figure  25  shows  the  isochromatics  calculated  from 
model  2  at  a  depth  of  10  km.  The  general  features  of  figure 
21  are  present  in  figure  25. 

Figure  26  shows  the  isopachs  calculated  from  model  2 
at  a  depth  of  10  km.  The  area  of  decreased  pressure  at  the 
surface  (figure  22)  parallel  to  the  White  Wolf  and  the 
Garlock  faults  is  replaced  by  an  area  of  slight  pressure 
increase,  but  the  other  general  features  seen  at  the  surface 
(figure  22)  are  present  10  km  deep. 

Figures  27  and  28  show  the  isochromatics  and  isopachs, 
respectively,  calculated  from  model  2  at  a  depth  of  20  km. 
Both  figures  show  the  same  general  features  as  figures  25 
and  26. 

Discussions  and  Conclusions 

Several  points  have  been  demonstrated  in  this  chapter. 
First,  a  north-south  oriented  compressive  stress  can  be 
produced  over  most  of  the  study  area  by  the  bend  in  the 
San  Andreas  fault.  This  is  important  because  it  removes 
the  necessity  of  invoking  such  a  regional  stress  system  to 
account  for  the  Transverse  Ranges  deformation.  Instead, 
north-south  compression  exists  as  a  natural  consequence 
of  the  bend  in  the  San  Andreas  fault. 

Second,  most  of  the  faults  in  the  Transverse  Ranges  that 
exhibit  Holocene  or  Quaternary  movement  are  parallel  or 
sub-parallel  to  the  calculated  slip  planes.  Thus,  the  stress 
field  produced  by  the  bend  in  the  San  Andreas  fault  proba- 
bly has  been  responsible  for  "selecting"  which  of  the  many 
faults  in  the  Transverse  Ranges  have  been  active.  For 
example,  by  reversing  the  movement  on  the  San  Andreas 
fault  by  50  to  75  km  (1  to  1.5  million  years  at  a  rate  of 
5  cm/yr),  the  western  parts  of  the  Santa  Ynez  and  the 
Arroyo  Parida  faults  would  be  moved  to  the  south  and 
east,  and  the  observed  movement  on  these  faults  would 
agree  with  that  predicted  by  model  2.  Recall  that  these 
faults  were  more  active  in  the  Pleistocene  than  in  the 
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Holocene  ("Structural  Characteristics";  Dibblee,  1966); 
this  is  to  be  expected  from  this  model. 

Third,  when  the  predicted  fault  patterns  of  models  1 
and  2  are  compared  to  the  observed  faults,  it  is  apparent 
that  model  2  predicts  the  observed  fault  patterns  better 
than  model  1.  This  rather  surprising  result  means  that  the 
depth  of  the  locked  zone  is  greater  than  would  have  been 
expected  based  on  data  from  central  California  (see 
"Mathematical  Model  of  the  San  Andreas  Fault"  and 
Appendix  A).  This  in  turn  indicates  that  past  geodetic 
measurements  (particularly  the  Geodimeter  measure- 
ments) may  have  been  made  too  close  to  the  San  Andreas 
fault  to  observe  strain  accumulation.  To  test  this,  the  Cali- 
fornia Division  of  Mines  and  Geology  intends  to  establish 
one  and  perhaps  two  trilateration  arcs  between  Tejon  and 
Cajon  Passes  that  will  extend  75  to  100  kms  approximate- 
ly normal  to  the  San  Andreas  fault. 

SEISMICITY  AND 

SECONDARY  FAULTING 

DETERMINED  FROM  THE  MODELS 

I  want  to  use  the  calculated  stresses  to  try  to  draw 
conclusions  about  (a)  the  occurrence  of  larger  earth- 
quakes in  the  study  area;  (b)  the  observed  faults  in  the 
study  area  which  have  evidence  of  Holocene  or  Quater- 
nary movement;  and  (c)  possible  zones  of  future  signifi- 
cant seismicity  in  the  study  area. 

Seismicity 
Areas  of  Potential  Failure 

Areas  of  potential  activity  can  be  determined  from  the 
isochromatic  and  isopach  maps  presented  earlier  if  failure 
criteria  are  established.  Chinnery  (1966a,  b)  discusses  a 
way  to  do  this  for  a  straight  strike-slip  fault.  He  first 
determines  the  shear  stress  drop  r  °,2  at  the  center  of  the 
fault;  then,  assuming  the  shear  stress  drop  relieves  all  the 
accumulated  shear  stress,  he  adds  two  different  regional 
stress  fields,  each  of  whose  t",,  component  equals  t  °,2  , 
to  the  calculated  stress  change  field.  The  resulting  stress 
fields  are  thus  the  stress  left  in  the  ground  after  an  earth- 
quake. He  then  assumed  that  ( 1 )  secondary  faulting  is 
likely  where  the  maximum  shear  stress  equals  or  exceeds 
T  °i2;  (2)  secondary  faulting  is  likely  in  regions  where  the 
average  compressive  stress  is  reduced;  and  (3)  secondary 
faulting  is  less  likely  when  the  average  compressive  stress 
is  increased. 

The  problem  in  the  present  work  is  somewhat  different 
in  that  I  have  assumed  there  is  no  pre-existing  regional 
stress  field.  Thus,  I  can  go  directly  to  the  calculation  of 
areas  of  potential  failure.  The  level  of  maximum  shear 
stress  necessary  to  produce  failure  must  first  be  deter- 
mined. Chinnery  (1966a)  determined  a  value  of  r  °,2  of 
about  5  bars  from  the  long  straight  strike-slip  fault  model 
be  used.  In  the  Transverse  Ranges,  faulting  is  more  often 
dip-slip  than  strike-slip,  so  this  value  of  t  °,;  might  not 


apply.  Thatcher  and  Hanks  (1973)  have  shown  that  stress 
drops  determined  from  S-wave  spectra  for  earthquakes  in 
the  Transverse  Ranges  with  magnitudes  between  2.9  and 
5.3  run  from  less  than  1  bar  to  about  80  bars,  with  essen- 
tially no  correlation  between  magnitude  and  stress  drop. 
Savage  and  Wood  ( 1971)  quote  a  stress  drop  of  about  10 
bars  for  the  magnitude  7.7  Arvin-Tehachapi  earthquake; 
and  estimates  of  the  stress  drop  of  the  magnitude  6.5  San 
Fernando  earthquake  run  from  21  to  34  bars  (Hanks, 
1974).  If  it  is  assumed  that  the  stress  drop  in  an  earth- 
quake completely  relieves  the  maximum  shear  stress 
(Tmaxl).  then  all  those  areas  where  r  max  exceeds  about 
1  bar  should  be  considered. 

From  the  above  discussion,  areas  of  potential  failure 
will  be  determined  by  the  following  rules: 

(a)  Foilure  is  most  likely  in  areas  where  the  average  compressive  stress  is 
reduced  and  where  the  maximum  shear  stress  is  greater  than  3  bors. 

(b)  Failure  is  likely  in  areas  where  the  average  compressive  stress  has  been 
slightly  increased  (by  about  1  bar)  and  the  maximum  shear  stress 
exceeds  5  bars. 

(c)  Failure  is  also  likely  in  areas  where  the  maximum  shear  stress  exceeds 
13  bars  and  the  increase  in  average  compressive  stress  does  not  exceed 
7  bars. 

For  rule  "a",  the  value  of  3  bars  is  between  the  1  bar 
or  less  found  by  Thatcher  and  Hanks  (1974)  and  the  5 
bars  used  by  Chinnery  (1966a).  If  the  average  compres- 
sive stress  has  been  increased,  it  seems  reasonable  to  re- 
quire failure  at  a  higher  value  of  the  maximum  shear 
stress.  Thus,  while  the  values  in  rule  "b"  are  somewhat 
arbitrary,  they  seem  reasonable.  Finally,  it  also  seems  rea- 
sonable to  say  that  failure  is  also  likely  to  occur  even  in 
those  regions  where  the  average  compressive  stress  has 
increased  if  the  maximum  shear  stress  is  very  high.  The 
values  of  13  bars  for  the  maximum  shear  stress  and  7  bars 
for  the  average  compressive  stress  are  also  arbitrary,  but 
they  do  not  appear  to  be  unacceptable. 

Historic  Seismicity  Compared  to 
Areas  of  Potential  Failure 

Using  these  rules  and  the  data  from  "Calculated  Stress 
Fields",  areas  of  potential  failure  were  determined  at  the 
surface  and  at  depths  of  10  and  20  km  for  both  models. 
Figures  29  to  32  show  these  areas  at  the  surface  and  at  a 
depth  of  20  km  for  models  1  and  2,  respectively.  The 
potential  failure  areas  at  the  10  and  20  km  depths  were 
similar  for  each  model,  and  hence,  only  the  potential  fail- 
ure areas  at  20  km  are  presented.  Figures  29  to  32  also 
show  the  epicenters  from  figure  5.  The  accuracy  of  the 
pre- 1932  epicenters  is  probably  no  better  than  ±  50  km, 
while  the  accuracy  of  the  post- 1932  epicenters  is  at  worst 
about  ±   15  km. 

The  only  events  that  do  not  lie  in,  or  near,  the  areas  of 
potential  failure  are  those  to  the  south  of  the  Malibu- 
Santa  Monica-Raymond  trend  and  three  events  offshore 
from  the  Arroyo  Parida  fault.  The  1910  and  1938  (and 
possibly  the  1903)  events  are  probably  associated  with  the 
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Elsinore  fault,  and  three  other  events  (1933  and  two  in 
1941)  are  assixiated  with  the  Newport-Inglewood  fault. 
Recall  that  both  the  Newport-Inglewood  and  the  Elsinore 
fault  zones  are  characterized  by  predominantly  right-lat- 
eral strike-slip  motion,  and  they  may  be  reacting  to  the 
same  stress  system  that  produces  movement  on  the  San 
Andreas  fault.  The  occurrence  of  the  other  events  in  the 
Santa  Barbara  Channel,  well  away  from  the  areas  of  po- 
tential failure,  is  more  difficult  to  account  for;  but  the 
events  are  probably  far  enough  from  the  San  Andreas  fault 
that  we  should  not  expect  the  models  to  have  much  effect. 

Note  that  the  events  in  the  eastern  comer  of  the  area  are 
generally  in  an  area  of  decreased  average  compressive 
stress,  but  the  maximum  shear  stress  is  less  than  3  bars. 
This  suggests  that  rule  "a"  in  "Areas  of  Potential  Failure" 
may  be  somewhat  conservative. 

Figures  29  to  32  suggest  that  the  areas  of  potential 
failure  derived  from  models  1  and  2  do  in  fact  delineate 
areas  of  high  seismicity  and  that  the  seismicity  in  the 
Transverse  Ranges  is  lower  outside  these  areas  of  potential 
failure.  Figures  31  and  32  further  suggest  that  model  2 
may  be  somewhat  better  than  model  1,  for  the  areas  of 
potential  failure  determined  from  model  2  include  some- 
what more  events  north  of  the  San  Andreas  fault. 

Secondary  Faulting 

Theory  suggests  that,  in  unfractured  material,  failure 
should  occur  on  the  slip  planes  (Jaeger  and  Cook,  1969). 
The  slip  planes  are  at  angles  of  ±  45°  to  the  largest  com- 
pressive stress  axis,  but  laboratory  evidence  indicates  that 
failure  of  initially  unfractured  rocks  usually  occurs  on 
planes  at  angles  of  between  30°  and  45°  to  the  largest 
compressive  stress  axis  (Jaeger  and  Cook,  1969).  If  the 
material  contains  fractures,  then  further  slippage  tends  to 
occur  on  the  previously  existing  fractures,  even  though 
these  fractures  may  not  make  angles  of  30°  to  45°  with  the 
axis  of  greatest  compression.  Thus,  the  slip  planes  will  be 
presented  as  the  simplest  way  of  visualizing  the  orienta- 
tion of  potential  slippage. 

By  combining  the  areas  of  likely  failure  (figures  29  to 
32)  with  the  slip  planes  presented  in  "Calculated  Stress 
Fields",  the  orientation  of  likely  fault  planes  can  be  deter- 
mined. Figure  33  shows  the  likely  slip  directions  deter- 
mined from  model  1  at  the  surface.  Of  the  left-lateral 
faults  showing  Holocene  and  Quaternary  movement  that 
lie  within  the  region  of  potential  failure,  only  the  Big  Pine 
fault  is  sub-parallel  to  a  likely  slip  direction.  A  few  right- 
lateral  faults  in  the  northwest  comer  of  figure  33  are 
within  the  region  of  potential  failure  and  have  trends  sub- 
parallel  to  likely  slip  directions.  The  extreme  northern  end 
of  the  San  Gabriel  fault  also  lies  sub-parallel  to  a  likely 
slip  direction. 

Figure  33  suggests  that  model  1  could  produce  the  ob- 
served sense  of  movement  on  a  few  of  the  faults  in  the 
study  area  in  the  upper  1  to  2  km  of  the  earth. 


Figure  34  shows  the  orientation  of  likely  fault  planes  at 
the  surface  determined  from  model  2.  The  left-lateral 
planes  are  sub-parallel  to  the  White  Wolf  and  the  Garlock 
faults  and  parts  of  the  Big  Pine  and  the  Santa  Ynez  faults. 
The  right-lateral  planes  are  sub-parallel  to  the  northem 
part  of  the  San  Gabriel  fault  and  to  the  faults  in  the 
northwestern  comer  of  the  figure. 

When  figures  33,  34,  15,  and  23  are  considered  together, 
it  is  apparent  that  model  2  more  accurately  predicts  the 
observed  sense  of  motion  on  the  White  Wolf  and  the 
Garlock  faults  than  does  model  1.  Model  2  also  delineates 
areas  of  higher  seismicity  more  clearly  than  does  model  1. 
While  neither  model  includes  the  1971  San  Fernando  epi- 
center within  an  area  of  potential  failure,  the  epicenter  of 
this  earthquake  does  lie  closer  to  the  area  of  potential 
failure  determined  from  model  2  than  that  determined 
from  model  1.  Thus,  the  results  of  this  study  suggest  that 
the  depth  of  the  locked  section  is  greater  than  would  have 
been  estimated  from  epicenter  data  ("Mathematical  Mod- 
el of  the  San  Andreas  Fault"). 

Discussion  and  Conclusions 

In  this  chapter,  we  have  shown  that  both  models  1  and 
2  define  reasonably  well  the  areas  in  the  Transverse 
Ranges  where  larger  earthquakes  have  occurred,  but  mod- 
el 2  tends  to  predict  the  observed  motion  on  mapped  faults 
better  than  model  1.  Thus,  I  can  use  these  results  with  a 
fair  degree  of  confidence  to  discuss  the  future  seismicity 
of  the  Transverse  Ranges. 

The  following  facts  must  be  kept  in  mind,  however.  The 
occurrence  of  any  earthquake  in  the  Transverse  Ranges 
does  alter  the  stress  fields  presented  in  this  Chapter;  and 
the  larger  the  earthquake,  the  greater  will  be  its  effect  on 
the  stress  field.  Thus,  to  arrive  at  the  "real"  stress  field  in 
the  Transverse  Ranges  at  the  present  time,  I  should  in- 
clude the  effects  of  all  the  earthquakes  which  have  oc- 
curred since  1857.  I  nonetheless  feel  that  the  results 
presented  here  can  be  used  to  draw  general  conclusions 
about  the  future  seismicity  of  the  Transverse  Ranges. 

Subject  to  these  reservations,  the  following  statements 
can  be  made.  Any  future  large  earthquakes  can  be  expect- 
ed to  occur  in,  or  very  near,  the  areas  of  potential  failure 
shown  on  figures  31  and  32.  Only  those  faults  that  trend 
parallel  or  sub-parallel  to  those  shown  will  be  involved  in 
the  future  large  earthquakes.  Of  course,  once  the  fault 
movement  has  been  initiated,  it  could  propagate  well 
beyond  the  areas  shown  in  figures  31  and  32. 

Specifically,  based  on  model  2,  future  large  earthquakes 
in  the  Transverse  Ranges  are  most  likely  on  the  Garlock 
fault,  the  Pleito  fault,  the  Frazier  fault,  the  Big  Pine  fault, 
and  the  northern  end  of  the  San  Gabriel  fault.  If  the  thrust 
faults  of  the  area — such  as  the  San  Cayetano,  the  Santa 
Susana,  and  the  Sierra  Madre — <:xtend  into  the  area  of 
potential  failure  at  depth,  large  earthquakes  can  be  expect- 
ed on  them.  Future  movement  on  the  Santa  Ynez,  the 
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Arroyo  Parida,  the  Malibu,  the  Santa  Monica,  and  the 
Raymond  faults  is  not  indicated  by  the  model.  Thus,  some 
stress  system  in  addition  to  that  produced  by  the  San 
Andreas  fault  must  be  invoked  to  account  for  the  observed 
Quaternary  and  Holocene  movement  on  these  faults. 

A  major  question  that  arises  is  whether  or  not  the  cal- 
culated stress  fields  are  large  enough  to  produce  the  ob- 
served earthquakes.  I  am  not  requiring  that  the  calculated 
stress  be  large  enough  to  cause  failure  of  virgin  rock;  but 
rather,  I  am  asking  that  the  stresses  be  large  enough  to 
initiate  sliding  on  previously  existing  fractures.  Intuitive- 
ly, I  would  expect  that  the  stress  required  to  initiate  sliding 
could  be  much  less  than  that  required  to  form  new  frac- 
tures (Jaeger  and  Cook,  1969).  Recall  also  that  the  stress 
drops  ( 1  to  80  bars)  determined  from  S-wave  spectra  for 
earthquakes  in  the  Transverse  Ranges  by  Thatcher  and 
Hanks  (1973)  and  by  Savage  and  Wood  (1971)  are  with- 
in the  range  of  calculated  values  of  maximum  shear  stress. 
The  deviatoric  horizontal  compressive  stresses  determined 
from  the  model  are  at  least  an  order  of  magnitude  less  than 
the  value  of  about  1  kilobar  determined  by  Hanks  (1974) 
for  the  1971  San  Fernando  earthquake.  However,  as 
Hanks  (1974)  points  out,  this  value  seems  to  apply  only 
to  the  first  seconds  of  the  fracture  process,  while  the  over- 
all stress  drop  is  on  the  order  of  30  bars.  It  is  entirely 
possible  that  this  high  stress  is  a  stress  concentration  due 
to  a  local  change  in  material  properties. 

It  also  should  be  pointed  out  that  this  model  cannot  be 
used  to  make  specific  earthquake  predictions.  The  model 
is  based  on  an  isotropic,  homogeneous  half-space  and  the 
Transverse  Ranges  are  anything  but  that.  Local  stress 
concentrations  undoubtedly  exist  wherever  the  material 
properties  change,  and  it  is  these  local  stress  concentra- 
tions that  must  be  responsible  for  individual  earthquakes. 

CONCLUSIONS  AND  SUGGESTIONS 
FOR  FURTHER  WORK 

The  vertical  displacement  and  stress  change  fields  have 
been  calculated  from  a  relatively  simple  dislocation  model 
of  the  San  Andreas  fault  in  the  Transverse  Ranges.  The 
calculated  vertical  displacement  fields  from  two  models 
resemble  the  general  topographic  features  of  the  Trans- 
verse Ranges,  and  this  suggests  that  the  present  topogra- 
phy may  have  been  produced  by  the  bend  in  the  San 
Andreas  fault.  It  does  not  appear  to  be  possible  to  select 
between  the  two  models  of  the  San  Andreas  fault  based  on 
the  topography  above. 

Several  conclusions  result  from  the  interpretation  of  the 
stress  change  fields  calculated  from  the  two  models. 

First,  it  does  appear  to  be  possible  to  select  between  the 
two  models  based  on  the  calculated  stress  fields.  Model  2, 
the  model  with  the  locked  zone  30  km  deep,  seems  to 
explain  the  observed  fault  patterns  and  strain  data  better 
than  model  1,  the  model  with  a  15-km-deep  locked  sec- 


tion. This,  of  course,  means  that  the  zone  of  stable  sliding 
is  deeper  in  the  Transverse  Ranges  than  in  central  Califor- 
nia (see  "Mathematical  Model  of  the  San  Andreas  Fault" 
and  Appendix  A). 

Second,  the  principal  compressive  stresses  calculated 
from  the  model  have  a  north-south  orientation  through- 
out most  of  the  study  area.  This  agrees  with  the  conclu- 
sions drawn  by  Hill  (1955),  Benioff  (1955),  and  Hanks 
(1974). 

Third,  the  areas  of  potential  failure  determined  from  the 
model  include  many  of  the  larger  earthquakes  in  the 
Transverse  Ranges;  and  the  larger  events  that  have  not 
occurred  in,  or  adjacent  to,  the  area  of  potential  failure  can 
be  related  to  movement  on  faults  in  the  Peninsular 
Ranges. 

Fourth,  the  slip  planes  determined  from  the  model  for 
the  1952  Arvin-Tehachapi  and  the  1971  San  Fernando 
earthquakes  have  azimuths  similar  to  those  of  the  ob- 
served fault  planes,  and  the  calculated  principal  stress 
system  could  certainly  produce  the  sense  of  movement 
observed  for  these  two  earthquakes. 

Finally,  the  calculated  stress  system  could  explain  the 
pattern  of  faults  having  Holocene  or  Quaternary  move- 
ment in  the  Transverse  Ranges.  This  does  not  necessarily 
imply  that  the  bend  in  the  San  Andreas  fault  is  responsible 
for  the  origin  of  the  Transverse  Ranges.  It  does,  however, 
mean  that,  given  a  complex  set  of  faults  of  many  different 
orientations,  the  stress  system  produced  by  the  bend  in  the 
San  Andreas  fault  will  cause  activity  mainly  on  those 
faults  that  are  favorably  oriented  to  the  stress  system  at 
the  time  of  observation.  For  the  same  time  interval,  other 
less  favorably  oriented  faults  will  show  a  much  lower  level 
of  activity,  even  though  these  faults  may  have  been  much 
more  active  at  some  earlier  time. 

A  complete  discussion  of  the  origin  of  the  Transverse 
Ranges  is  beyond  the  scope  of  this  paper.  However,  some 
brief  comments  must  be  made.  The  models  discussed  here 
do  not  produce  stress  fields  that  are  large  enough  to  initi- 
ate failure  in  fresh  rock.  Thus,  the  bend  in  the  San  An- 
dreas fault  probably  did  not  produce  the  Transverse 
Ranges.  But  I  do  feel  that  the  stresses  produced  by  the 
bend  are  large  enough  to  niaintain  at  least  some  of  the 
Transverse  Ranges  structures  after  these  structures  enter 
the  stress  fields. 

There  is  reasonably  good  evidence  for  much  greater 
displacements  on  the  San  Andreas  fault  system  north  of 
the  Transverse  Ranges  than  on  the  fault  system  south  of 
the  Transverse  Ranges  (Berry,  1973;  Dickinson  and  oth- 
ers, 1972).  The  interpretation  of  this  evidence  is  that  the 
San  Andreas  fault  system  south  of  the  Transverse  Ranges 
is  younger  than  the  northern  part  of  the  system  and  that 
it  probably  took  some  time  for  the  fault  to  propagate  south 
from  the  present  location  of  the  Transverse  Ranges.  Thus, 
the  fault  probably  "ended"  in  the  vicinity  of  the  Tehach- 
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api  Mountains  for  a  considerable  length  of  time,  and  the 
analysis  used  by  Chinnery  (1966a,  b)  could  be  modified 
and  applied  to  attempt  to  account  for  some  or  all  of  the 
structures  of  the  Transverse  Ranges.  This  seems  to  be  a 
promising  line  of  attack  on  the  problem  of  the  origin  of  the 
Transverse  Ranges. 

It  is  rather  remarkable  that  the  simple  model  of  a  dislo- 
cation in  an  isotropic,  homogeneous,  linearly  elastic  half- 
space  can  provide  a  fairly  reasonable  approximation  of  the 
observed  faulting  in  such  a  heterogeneous  environment  as 
the  Transverse  Ranges.  The  model  can  be  improved  and 
made  more  realistic  in  several  ways.  First,  the  assumption 
of  an  elastic  half-space  can  be  relaxed  and  non-linear 
effects  included  by  using  the  work  of  Rosenman  and  Singh 
(1973)  and  Singh  and  Rosenman  (1974).  This  would  not 
involve  major  changes  to  the  computer  programs  used 
here,  and  it  would  allow  a  more  realistic  time  history  of 
the  stress  accumulation  to  be  studied. 

At  the  same  time,  the  results  of  Mansinha  and  Smylie 
( 1971 )  can  be  incorporated  to  study  the  effects  that  move- 
ment on  some  of  the  thrust  faults  in  the  area  would  have 
on  the  stress  system  due  to  the  bend  in  the  San  Andreas 
fault.  Thus,  we  would  be  able  to  follow  the  stress  field 


through  several  cycles  of  accumulation  and  release  on  the 
San  Andreas  fault  and  perhaps  study  the  development  of 
the  Transverse  Ranges.  These  studies  would  still  be  an 
isotropic,  homogeneous  half-space  model  and  would  be 
fairly  simple  to  carry  out,  for  the  equations  all  exist  in  a 
closed  form  and,  hence,  can  be  readily  evaluated. 

To  construct  an  even  more  realistic  model,  a  three- 
dimensional  finite  element  model  of  the  region  would  have 
to  be  used.  This  would  mean,  of  course,  that  the  basic 
equations  would  have  to  be  integrated  numerically  for 
each  model  one  is  interested  in.  This  would  very  likely  be 
a  time-consuming,  expensive  task,  even  on  the  present 
generation  of  high-speed  computers. 

Finally,  any  further  work  with  these  models  will  be 
hampered  by  the  lack  of  geodetic  observations  in  the  re- 
gion. Thus,  it  will  be  very  difficult  to  select  one  model  as 
"better"  than  another.  As  a  result  of  the  present  study,  the 
California  Division  of  Mines  and  Geology  is  redirecting 
part  of  its  Geodimeter  measurement  program  to  be  able 
to  obtain  a  greater  amount  of  relevant  data.  In  particular, 
trilateration  arcs  extending  50  to  75  km  from  the  fault  will 
be  observed  in  an  attempt  to  obtain  more  data  on  fault 
movements  at  depth. 
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Appendix  A 

STRAIN  TRANSITION  ZONE  ON  THE  SAN  ANDREAS  FAULT 
BETWEEN  PARKFIELD  AND  TAFT,  CALIFORNIA 


Several  authors  have  pointed  out  that  fault  creep  on  the  San 
Andreas  fault  dies  out  in  the  vicinity  of  Parkfield  (Allen,  1968; 
Nason,  1971;  Yamashita  and  Burford,  1973),  and  the  rapid 
decrease  in  seismicity  south  of  Parkfield  has  also  been  noted 
(Brune  and  Allen,  1967;  Allen,  1968;  "Seismicity",  this  work). 
Allen  (1968)  suggested  that  the  segment  of  the  San  Andreas 
fault  south  of  Parkfield  is  locked,  and  Rodgers  and  Chinnery 
(1973)  presented  a  model  of  the  San  Andreas  fault  incorporating 
the  locked  section. 

There  were  exhaustive  studies  of  the  continuing  deformation 
along  the  fault  for  several  months  following  the  1966  Parkfield 
earthquakes;  but  to  my  knowledge,  no  work  has  been  published 
that  considers  the  details  of  strain  variation  along  the  San  An- 
dreas fault  since  then.  In  this  work,  I  will  examine  in  some  detail 
the  geodetic  evidence  south  of  Parkfield  that  delineates  the  strain 
transition  zone,  and  I  will  present  a  dislocation  model  that  ade- 
quately explains  the  major  findings  of  the  geodetic  work. 

PREVIOUS  WORK 

Cherry  and  Savage  (1972)  analyzed  the  results  of  repeated 
trilaterations  of  the  pentagonal  survey  figure  near  Parkfield  (fig- 
ure A-1).  They  fitted  straight  lines  to  the  line-length  changes, 
plotted  the  rate  of  length  change  against  the  azimuth  of  the  line, 
and  fitted  a  strain  field  to  the  results.  They  found  a  strain  field 
with  an  essentially  east-west-oriented  principal  extension  and 
very  little  or  no  compression.  Greensfelder  and  Bennett  (1973) 
suggested  that  a  straight  line  fit  to  the  data  does  not  reflect  what 
is  actually  occurring  and  that  a  sinusoidal  curve  is  a  better  fit  to 
the  data. 

GEODETIC  EVIDENCE  OF  THE  TRANSITION  ZONE 

Figure  A-1  shows  selected  lines  from  the  California  Geodime- 
ter  Network  and  all  operating  creepmeters  south  of  Stone  Can- 
yon (XSCl).  It  also  shows  the  location  of  fault-crossing  survey 
site  TEM  (see  Appendix  B).  The  area  I  shall  be  most  concerned 
with  is  between  Stone  Canyon  and  Simmler. 

Let  me  first  consider  the  creepmeter  data  (table  A-1).  The 
creepmeter  at  Carr  Ranch  may  not  be  across  the  presently  active 
fault  trace,  which  would  explain  the  zero  creep  rate  between 
XDR 1  and  XGH I .  It  is  clear  that  the  creep  changes  dramatical- 
ly between  XDRl  and  Twisselman  Ranch,  going  to  essentially 
zero  over  a  distance  of  about  30  km.  Again,  it  is  possible  that  the 
instruments  south  of  XDRl  do  not  cross  the  presently  active 
fault  trace;  but,  since  the  U.S.  Geological  Survey  instruments 
were  installed  across  the  surface  break  of  the  1966  Parkfield 
earthquakes  (Yamashita  and  Burford,  1973),  this  does  not  seem 
very  likely.  Also,  creep  is  observed  at  XGHl  and  XWTl;  and, 
even  though  the  rate  is  very  low,  this  indicates  that  the  fault  trace 
on  which  the  instruments  are  located  is  at  least  intermittently 
active. 


Table  A-1.  Creepmeter  data.  All  instruments  installed  in  1969  and 
still  operating  in  1974. 


Site 


Creep 

rate 

(mm/yr) 


Source 


XSCl 

24 

XDRl 

11 

Carr  Ranch 

0 

XGHl  2 

XWTl  3 

Twisselman  Ranch    0 


Yamashita  and  Burford  (1973) 
Yamashita  and  Burford  (1973) 
James  Hileman,  California  Institute  of 

Technology,  personal  communication, 

1974 
Yamashita  and  Burford  (1973) 
Yamashita  and  Burford  (1973) 
James  Hileman,  California  Institute  of 

Technology,  personal  communication, 

1974 


The  survey  figure  TEM  (figure  A-1)  is  discussed  in  some 
detail  in  Appendix  B.  There  is  a  suggestion  of  right-lateral  dis- 
placement across  the  fault  with  a  rate  of  several  millimeters  per 
year;  but,  as  discussed  in  Appendix  B,  this  may  be  an  effect  due 
to  the  rapid  change  in  the  creep  rate  in  this  area. 

Observations  of  the  California  Goedimeter  Network  also  give 
some  indication  of  the  extent  of  the  transition  zone  (figures  A-2, 
A-3,  and  A-4).  The  dashed  lines  on  the  plots  are  included  only 
to  suggest  trends  in  the  data.  They  are  fitted  to  the  data  by  eye. 

Lines  25  through  31  in  figure  A-2  show  a  reasonably  smooth 
change  in  length  from  1959  through  1973.  The  average  rate  for 
these  lines  is  about  27  mm /year  (Savage  and  others,  1973). 
Note  the  abrupt  change  in  rate  in  line  32  before  and  after  the 
1966  Parkfield  earthquakes.  The  reason  for  this  change  is  not 
clear,  for  line  31,  with  nearly  the  same  azimuth,  does  not  show 
such  a  change  in  rate.  However,  the  important  point  here  is  the 
essentially  linear  trends  shown  by  all  the  lines  north  of  Parkfield. 

Figure  A-3  shows  Geodimeter  lines  between  Parkfield  and 
Maricopa  which  cross  the  fault  at  acute  angles.  Note  that  line 
33  has  a  linear  trend  from  1959  to  1963  and  from  1966  to  1969. 
It  has  been  suggested  (Hofmann,  1968;  Greensfelder  and  Ben- 
nett, 1973)  that  the  anomalous  change  in  the  line  length  between 
1964  and  1966  was  a  precursor  to  the  1966  earthquakes.  If  this 
is  true,  then  the  line  behavior  after  1969  suggests  that  an  earth- 
quake may  be  imminent  in  the  Parkfield  area. 

Prior  to  1966,  lines  35  and  39  show  very  little  net  change  in 
line  length.  After  the  1966  Parkfield  earthquakes,  both  lines 
show  a  net  change  in  length  of  several  centimeters;  and,  from 
about  1969  on,  the  two  lines  show  very  little  net  change.  Lines 
40  and  46  (figure  A-3)  show  essentially  no  change  from  1959 
to  1973. 
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Thus,  the  Geodimeter  hnes  that  cross  the  fault  at  an  angle 
suggest  that  there  is  a  change  in  the  strain  pattern  on  the  segment 
of  the  fault  spanned  by  hnes  33  and  35. 

Similarly,  Geodimeter  lines  approximately  perpendicular  to 
the  fault  (figure  A^)  indicate  a  change  in  the  strain  pattern 
between  lines  291  and  297.  Measurements  were  initiated  on  line 
297  to  fill  in  between  lines  291  and  263.  Unfortunately,  soon  after 
line  297  was  established,  the  large  changes  in  line  263  started 
dying  out.  Note  that  the  behavior  of  line  297  does  agree  with  that 
of  line  263. 

To  summarize,  the  geodetic  data  indicate  a  major  change  in 
the  strain  release  pattern  along  the  San  Andreas  fault  between 
Parkfield  and  Simmler  (figure  A-I).  The  data  from  measure- 
ments very  near,  or  on,  the  fault  trace  suggest  that  this  change 
occurs  between  TEM  and  the  creepmeter  at  Twisselman  Ranch. 


Using  the  model  shown  in  figure  A-5,  the  changes  in  Geodim- 
eter line  lengths,  dl/dt,  were  calculated  for  the  lines  shown  in 
figures  A-5  and  A-6.  The  straight  lines  on  the  graphs  show  the 
calculated  length  change  rates. 

In  addition,  a  displacement  perpendicular  to  the  fault  of  the 
form 

U  =  U„  sin  (t/T) 

with  Uo  =  3  cm,  t  =  1  years,  and  t  in  years  was  added  to  the 
displacement  field  calculated  from  the  model.  This  is  shown  by 
the  sinusoidal  curves  in  figures  A-5  and  A-6.  The  phase  was 
determined  such  that  the  best  fit  was  obtained  for  line  F  in  figure 
A-6.  Thus,  all  the  changes  are  in  phase.  The  period  of  seven 
years  is  the  average  for  the  whole  network  as  determined  by 
Greensfelder  and  Bennett  (1973).  The  amplitude  of  3  cm  is 
somewhat  arbitrarily  chosen. 


DISLOCATION  MODEL  OF  THE  TRANSITION  ZONE 

The  dislocation  model  we  will  construct  must  satisfy  the  geo- 
detic observations  of  the  previous  section.  In  addition,  the  model 
will  try  to  account  for  the  periodic  changes  in  Geodimeter  line 
length  discussed  by  Greensfelder  and  Bennett  (1973). 

The  displacement  fields  calculated  from  several  possible  mod- 
els were  compared  to  the  Geodimeter  observations.  The  model 
which  best  fit  the  data  is  fairly  simple  and  will  be  presented  here. 
The  model  starts  with  an  infinitely  long  straight  fault  slipping 
from  20  km  to  infinity  at  a  rate  of  3  cm/year  (figure  A-5).  The 
displacement  field  is  calculated  using  the  expressions  given  by 
Savage  and  Burford  (1973). 

To  this  is  added  the  displacement  field  of  the  fault  surface 
shown  in  figure  A-5  which  has  a  0.5-km-deep  locked  segment 
from  Parkfield  to  near  line  291,  and  a  20-km-deep  locked  seg- 
ment from  line  291  south.  A  model  was  tried  that  had  a  4-km- 
long,  0.5-km-deep  slipping  segment  centered  on  TEM,  but  it  did 
not  significantly  improve  the  fit  of  the  model  to  the  data. 


DISCUSSION  AND  CONCLUSIONS 

The  calculated  length  change  rates  in  figure  A-6  agree  well 
with  the  observed  data  except  for  the  strain  steps  due  to  the  1966 
Parkfield  earthquakes.  The  sinusoidal  curves  do  not  significantly 
improve  the  fit  of  the  calculated  length  change  rates  to  the  data 
for  the  long  Geodimeter  lines. 

The  linear  rates  calculated  for  the  Parkfield  pentagon  do  not 
fit  the  data  very  well;  but  except  for  line  A,  the  sinusoidal  rates 
give  a  remarkably  good  fit  to  the  data  (figure  A-6).  And  even 
for  line  A,  the  calculated  length  changes  agree  in  phase  but  not 
in  slope. 

Thus,  we  have  derived  a  simple  dislocation  model  that  ac- 
counts for  the  trends  in  the  Geodimeter  data.  In  addition,  we 
have  indicated  that  the  periodic  fluctuations  in  line  length  point- 
ed out  by  Greensfelder  and  Bennett  ( 1973)  can  be  accounted  for 
by  a  simple  assumption.  The  origin  of  a  periodic  displacement 
perpendicular  to  the  San  Andreas  fault  is  not  clear  at  this  time. 
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EXPLANATION  OF  SYMBOLS 

^      Geodlmeter  Stotlons 

O  Towns 

O  Creepmeters 

'-'  Survey  Figure 

19]  Geodlmeter  Line  Number 
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Figure  A-1.        Creepmeters,  Geodlmeter  lines,  and  small  survey  figure  in  the  vicinity  of  Parkfield. 
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Figure  A-3.        Geodimeter  lines  and  data,  1959-73, 
Parkfield  to  Maricopa. 
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Figure  A~4.        Geodimeter  lines  normal  to  the  San 
Andreoi  fault,  1971-73,  Parkfield  to  Maricopo. 
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Figure  A-5.        Fault  model  and  predicted  and  observed  line  length  changes. 
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Figure   A-6.        Predicted   and   observed   line   length  changes,  Parkfield  pentagon. 


1979 


DEFORMATION.  STRESS  AND  FAULTING  TRANSVERSE  RANGES 


59 


Appendix  B 
FAULT-CROSSING  SURVEY  FIGURES 


INTRODUCTION 

In  1964-1966,  thirty  small  survey  figures  were  established 
across  several  faults  of  the  San  Andreas  fault  system  by  the  U.S. 
Coast  and  Geodetic  Survey  in  cooperation  with  the  California 
Department  of  Water  Resources  (Meade,  1966,  1968,  1971a, 
1971b).  The  data  and  preliminary  analyses  of  results  from  reob- 
servation  of  these  figures  have  been  released  in  an  informal 
manner  in  a  series  of  repwrts  titled  "Results  of  Triangulation  for 
Earth  Movement  Study  at  California  Aqueduct-Fault  Crossing 
Sites".  Meade  (1966,  1968,  1971a,  1971b)  briefly  reported  on  the 
results  from  these  survey  figures.  However,  to  my  knowledge,  a 
systematic  comparison  and  interpretation  of  these  results  has  not 
been  pubhshed.  Hence,  in  this  Appendix,  I  will  present  the  re- 
sults from  the  surveys  of  these  figures  and  attempt  a  partial 
analysis  of  the  results. 

Figure  B-1  shows  the  locations  of  the  survey  sites.  I  will 
discuss  each  site,  moving  from  north  to  south. 

TOLAY 

This  figure  was  established  on  the  north  shore  of  San  Pablo 
Bay  in  1965  and  resurveyed  once,  in  1967  (Greensfelder,  1972). 
The  figure  was  supposed  to  span  the  possible  northern  extension 
of  the  Hayward  fault,  but  the  exact  location  of  any  fault  in  this 
area  is  difficult  to  ascertain  (Rice,  1974).  In  addition,  there  were 
large,  inexplicable  triangle  closures  in  the  1965  survey  (Miller, 
1967a).  First-order  levels  were  also  run  at  the  site  in  1965. 

All  the  angle  changes  were  within  survey  error  and,  hence,  any 
fault  movement  is  unresolvable. 

MAR  (MIRA  VISTA) 

This  site  was  established  on  the  Hayward  fault  (figure  B-2) 
in  1966  and  resurveyed  in  1968  and  1970.  Miller  ( 1970c)  notes 
that  the  1968  and  1970  results  are  probably  more  reliable  since 
the  astronomic  azimuths  in  those  years  are  in  good  agreement, 
while  the  1966  azimuth  is  some  5  seconds  smaller. 

There  is  some  indication  of  right-lateral  movement  shown  in 
the  position  vectors,  but  Miller  (1970c)  suggests  this  may  be  due 
to  changes  in  the  astronomic  azimuths.  It  is  interesting  to  note 
that  the  astronomic  azimuths  of  both  B-A  and  B-C  increased. 
Tliis  suggests  there  could  be  another  unrecognized  fault  trace 
between  stations  B  and  C. 

STADIUM 

This  site  was  established  on  the  nm  of  Memorial  Stadium  on 
the  University  of  California's  Berkeley  campus  in  1966.  Earlier 
observations  had  indicated  that  the  stadium  was  built  across  the 
Hayward  fault  (figure  B-3)  and  that  creep  on  the  fault  was 
pulling  it  apart  (Bolt  and  Marion,  1966).  The  taped  base  line 


and  the  observed  astronomic  azimuth  of  the  original  survey  in 
1966  agree  with  those  observations  in  the  resurveys  of  1967  and 

1968  (Miller,  1970c).  Thus,  this  figure  can  be  considered  to  be 
quite  stable. 

Changes  in  azimuth  of  lines  A-D  and  B-C  suggest  an  average 
creep  rate  over  the  three  years  of  5  mm/year  (Miller,  1970c). 

CAMP  PARKS 

This  site  consists  of  four  figures  established  in  1964.  Three  of 
the  figures  cross  the  Pleasanton  fault  zone  (figure  B-4),  and  the 
fourth  is  a  line  across  the  Calaveras  fault.  I  will  not  discuss  the 
alignment.  The  figures  were  resurveyed  in  1965  and  1969  (Mil- 
ler, 1970d).  The  1969  survey  tied  the  three  figures  on  the  Pleas- 
anton fault  together.  The  maximum  triangle  closures  are  within 
first-order  survey  specifications.  However,  in  the  figure  PFS,  the 
base  line  2-5  increased  4  mm  and  underwent  an  astronomic 
azimuth  change  of  4  seconds  between  1964  and  1965.  The  length 
increase  of  the  base  line  would  account  for  the  overall  expansion 
of  the  figure  in  1965.  The  base  line  and  its  azimuth  remained 
relatively  constant  in  1969,  thus  suggesting  that  the  observed 
length  changes  may  be  more  reliable.  The  length  and  azimuth  of 
the  base  hne  3-4  of  PFE  suggest  that  the  observed  line-length 
changes  are  fairly  reliable.  The  base  line  1-2  of  PFW  underwent 
a  shortening  of  1  cm  between  1965  and  1969,  which  would 
account  for  the  large  amount  of  shortening  of  the  figure  between 
1965  and  1969. 

Miller  (1970d,  p.  4)  interprets  these  results  as  follows: 

"At  Site  PFS,  there  was  systematic  right-lateral  movement  of  about  18  mm 
for  the  5-year  interval  1964-1969.  Chonges  at  Site  PFE  do  not  indicote 
significant  movement,  but  at  Site  PFW  from  1965  to  1969,  there  is  on 
indication  of  about  1  cm  relative  movement." 

However,  the  changes  in  the  astronomic  azimuths  of  the  base 
lines  may  indicate  regional  strain  accumulation  rather  than  fault 
shppage.  The  three  figures  each  have  two  stations  in  the  fault 
zone,  and  this  may  introduce  unpredictable  changes  into  the 
results.  It  is  probably  necessary  in  a  situation  like  this  to  have 
several  longer  lines  that  would  span  the  whole  fault  zone  to 
resolve  the  problem. 

In  1974,  personnel  from  the  California  Division  of  Mines  and 
Geology  observed  these  figures,  using  an  MA-100  Tellurometer. 
The  observed  line  lengths  agreed  with  the  adjusted  length  of  the 

1969  triangulation  survey  to  within  the  errors  of  the  Tellurome- 
ter (±  1.5  mm).  The  more  recent  work  suggests  that,  if  fault 
slippage  occurred,  it  occurred  as  an  event  between  1965  and 
1969.  This  would  be  in  agreement  with  creepmeter  data  which 
shows  that  fault  creep  on  the  Hayward-Calaveras  fault  is  more 
"eventful"  than  steady-state  (Nason  and  others,  1974). 
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Figure  B-1.        Survey  figurei  ot  oqueducl-fault  crossing  sites   (modified  from  Greensfelder,  1972,  p.  8) 
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GREEN 

This  site  was  established  on  the  Greenville  fault  in  1965  and 
resurveyed  in  1966.  First -order  levels  were  run  at  the  site  in  1965 
and  1967.  The  astronomic  azimuth  of  the  line  C-A  remained 
essentially  constant,  but  the  taped  length  A-C  changed  by  1  cm. 
Miller  (1966b)  says  that  the  change  in  length  of  A-C  was  not 
supported  by  any  angle  changes  and  thus  is  probably  due  to  the 
difficulty  of  taping  this  line.  Triangle  closures  and  angle  correc- 
tions were  large,  but  they  were  within  the  specifications  of  the 
survey.  There  does  not  appear  to  be  any  consistent  horizontal  or 
vertical  movement  at  this  site. 

UNION 

This  site  on  the  Hay  ward  fault  (figure  B-5)  was  established 
in  1965  and  resurveyed  in  1966,  1967,  and  1969.  Triangle  clo- 
sures and  angle  corrections  were  within  first-order  specifica- 
tions. The  astronomic  azimuth  of  the  base  line,  A-C,  decreased 
by  less  than  one  second.  The  taped  length  A-C  changed  by  8  mm 
between  1966  and  1967.  First-order  levels  were  run  in  1965  and 
1967. 

Changes  in  directions  of  the  fault-crossing  lines  suggest  a 
creep  rate  of  about  8  mm/year  on  the  Hay  ward  fault  (Miller, 
1970a).  Changes  in  elevation  indicate  stations  west  of  the  fault 
rose  relative  to  stations  east  of  the  fault. 

VERAS 

This  site  was  established  in  1965  on  the  Calaveras  fault  (figure 
B-6)  and  resurveyed  in  1966  and  1970.  Miller  (1970e)  notes, 
"The  triangle  closures  for  the  surveys  of  1965  and  1966  are  larger 
than  normally  expected,  but  within  the  accuracy  limits  of  the 
specification..."  The  astronomic  azimuth  and  taped  length  of  the 
base  C-E  did  not  change  significantly  between  the  surveys. 
First-  order  levels  were  run  in  1965  and  1967. 

Again,  the  fault-crossing  lines  show  a  consistent  change  in 
direction  that  could  be  caused  by  right-lateral  creep  on  the 
Calaveras  fault.  Miller  (1970e)  says,  "Overall  observed  changes 
in  the  horizontal  directions  of  lines  approximately  perpendicular 
to  the  fault  showed  an  average  right-lateral  movement  of  12  mm 
for  the  5-year  interval."  The  elevation  differences  indicate  that 
marks  east  of  the  fault  rose  relative  to  marks  west  of  the  fault. 

OSGOOD 

This  figure  was  established  across  the  Hayward  fault  (figure 
B-7)  in  1966  and  resurveyed  in  1968  and  1969.  Triangle  closures 
and  angle  corrections  were  within  first-order  specification.  The 
taped  length  of  A-C  changed  by  3  mm,  but  the  astronomic 
azimuth  of  the  line  A-C  changed  nearly  4  seconds,  which  sug- 
gests that  there  may  be  some  strain  accumulation  in  the  vicinity 
of  the  net. 

Each  of  the  lines  crossing  the  fault  trace  shows  a  net  increase 
in  azimuth  over  four  years.  Resolving  these  angles  into  move- 
ment parallel  to  the  fault,  a  creep  rate  of  about  7  mm/year 
(Miller,  1970c)  is  obtained.  This  is  consistent  with  the  rate  of 
5  mm/year  at  STADIUM  and  8  mm/year  at  UNION. 

CIENEGA  VALLEY   (HARRIS,  TAYLOR,  STONE) 

These  figures  have  been  established  in  Cienega  Valley;  the 
earliest  was  The  resurveys  at  these  sites  indicate  a  creep  rate  of 
between  10  and  15  mm/year.  Meade  (1971  a,b)  has  discussed 
these  figures,  so  no  further  discussion  will  be  given  here. 


TEM 

This  figure  was  established  across  the  San  Andreas  fault  (fig- 
ure B-8)  in  1964  and  resurveyed  in  1965,  1966,  1968,  and  1971. 
All  the  triangle  closures  and  corrections  were  within  the  limits 
of  the  survey  specifications.  The  astronomic  azimuth  of  line  B-D 
remained  essentially  constant,  but  the  taped  length  of  B-D 
changed  by  12  mm  between  1966  and  1971  (Miller,  1972a). 
Levels  were  run  in  1968  and  1971. 

There  was  no  significant  difference  between  the  1964  and  1965 
surveys;  but  the  1966  survey,  made  after  the  1966  Parkfield 
earthquakes,  indicated  a  25  mm  right-lateral  shift  across  the 
fault.  The  1968  and  1971  surveys  suggest  continuing  right-lat- 
eral movement  at  a  rate  of  about  7  mm/year  (Miller,  1972a). 
The  elevation  differences  indicate  that  all  the  marks  are  sinking 
relative  to  mark  F,  but  the  three  westernmost  marks  (B,  C,  and 
D)  are  sinking  faster  than  the  other  marks. 

Detailed  mapping  following  the  1966  earthquakes  indicated 
that  ground  surface  rupture  ended  some  4.5  km  northwest  of 
TEM  (Brown  and  Vedder,  1967).  Eaton  and  others  (1970) 
show  that  the  aftershock  zone  of  the  1966  earthquakes  terminat- 
ed abruptly  about  5  km  northwest  of  TEM.  Thus,  the  ijf)served 
azimuth  changes  at  TEM  in  1966  could  reflect  either  movement 
on  a  fault  plane  that  did  not  intersect  the  surface  or  the  end 
effects  due  to  a  dislocation  surface  such  as  Chinnery  (1961, 
1963)  discussed.  The  fact  that  the  aftershock  zone  did  not  extend 
to  the  vicinity  of  TEM  suggests  that  the  azimuth  changes  are 
probably  due  to  end  effects  rather  than  movement  on  a  fault 
plane. 

The  two  resurveys  since  1966  suggest  continuing  fault  move- 
ment after  the  1966  earthquakes.  However,  data  from  the  U.S. 
Geological  Survey  microearthquake  net  does  not  show  any  seis- 
mic activity  south  of  the  vicinity  of  Parkfield.  The  rather  large 
change  (12  mm)  in  the  length  of  the  base  line  B-D  suggests 
strain  accumulation  rather  than  fault  creep.  Thus,  the  result  of 
the  TEM  resurveys  indicate  either  aseismic  fault  creep  that  prob- 
ably does  not  intersect  the  ground  surface  or  end  effects  due  to 
the  termination  of  surface  creep  between  creepmeters  XGHl 
and  XWTl  (Appendix  A).  Since  the  azimuth  changes  are  large 
even  for  lines  that  do  not  cross  the  fault,  neither  hypothesis  can 
be  ruled  out. 

SANTA 

This  site  was  established  on  the  Nacimiento  fault  (figure  B-9) 
in  1964  and  reobserved  in  1965  and  1968.  The  triangle  closures 
and  angle  corrections  were  large  in  all  three  surveys.  Astronomic 
azimuths  of  the  line  A-C  remained  essentially  constant,  and  the 
taped  length  of  A-C  lengthened  by  4  mm  between  1964  and 
1965.  The  survey  results  do  not  indicate  movement  across  the 
Nacimiento  fault  larger  than  the  survey  errors  (Miller,  1968). 
A  level  survey  was  run  in  1968,  and  the  levels  have  not  been 
rerun  since  then. 


METTLER 

This  site  was  established  in  1964  across  the  White  Wolf  fault 
and  resurveyed  in  1965,  1966,  and  1969.  There  is  good  evidence 
that  several  of  the  marks  in  the  figure  were  disturbed  between 
the  resurveys,  and  the  site  has  been  abandoned  (Miller.  1970b). 
Due  to  the  problems  with  the  site,  I  feel  that  the  results  of  the 
resurveys  are  of  such  poor  quality  that  no  further  analysis  is 
warranted. 
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MEADE 

This  site  was  established  (figure  B-10)  in  1969  to  replace 
METTLER  and  resurveyed  in  1972.  First-order  levels  were  run 
on  the  site  in  1969,  1970.  and  1972.  Triangle  closures  and  angle 
corrections  were  within  the  specifications  of  the  survey,  and  the 
non-fault-crossing  taped  lines  (A-B  and  D-F)  showed  no  sig- 
nificant changes.  However,  the  a.stronomic  azimuth  ofline  A-F 
decreased  by  3  seconds  (Miller,  1972d). 

The  adjusted  azimuth  changes  are  dominated  by  a  change  in 
the  astronomic  azimuth  ofline  A-F.  Miller  (1972d)  states  that 
the  location  of  the  White  Wolf  fault  is  not  well  known  at  this  site. 
The  counterclockwise  rotation  ofline  A-F  and  the  depression  of 
marks  A  and  B  relative  to  RON  would  be  consistent  with  move- 
ment on  the  White  Wolf  fault  if  it  passed  through  the  northern 
part  of  the  survey  figure. 

RANCH-TEJON 

In  1964,  RANCH  and  TEJON  were  established  on  the  Gar- 
lock  fault  (figures  B-1 1  and  B-12).  They  were  both  resurveyed 
in  1965  and  1966;  and  RANCH  was  resurveyed  in  1967,  1969, 
1970,  and  1972.  In  1967.  most  of  the  marks  at  the  TEJON  site 
were  destroyed  by  construction  on  the  California  Aqueduct.  In 
1970,  a  new  set  of  marks  was  established  near  the  original  TE- 
JON site,  and  RANCH  and  TEJON  were  connected  by  several 
longer  lines. 

RANCH 

Levels  were  run  at  RANCH  site  in  1964.  1965,  1966,  1968, 
1969.  1971.  and  1972;  and  table  B-1  shows  the  elevation  differ- 
ences. The  elevation  changes  indicate  that  all  marks  were  rising 
relative  to  D  from  1964  to  1968;  in  the  interval  1968  to  1969.  all 
the  marks  subsided  relative  to  D.  and  from  1969  to  1972.  the 
marks  generally  moved  up  relative  to  D.  but  the  amount  of 
change  was  less  than  in  earlier  periods.  The  results  of  the  hori- 
zontal surveys  show  that  all  lines  crossing  the  fault  except  C-F 
lengthened,  while  C-F  shortened  by  a  small  amount.  The  non- 
fault-crossing  lines  all  shortened  substantially,  reflecting  the 
shortening  of  the  taped  base  line  D-E.  Thus,  the  adjusted  line 
lengths  suggest  an  expansion  across  the  fault  with  some  left- 
lateral  movement.  The  observed  directions  from  the  survey  indi- 
cate left-lateral  movement  at  a  rate  of  between  6  and  8  mm  /year 
across  the  Oarlock  fault  (Miller.  1970b.  1972d).  If  this  is  a  real 
measure  of  fault  creep,  it  indicates  that  the  Oarlock  fault  near 
this  site  is  one  of  the  most  active  faults  in  California. 


However,  a  note  of  caution  should  be  sounded.  Between  1964 
and  1968.  a  tunnel  was  being  driven  between  RANCH  and 
TEJON.  and  the  south  portal  of  the  tunnel  is  in  the  middle  of 
the  RANCH  figure.  Thus,  the  horizontal  and  vertical  changes 
noted  above  might  be  a  result  of  the  tunneling.  However,  the 
consistency  of  the  observations  argues  against  there  being  a  large 
effect  at  the  site  due  to  the  tunnel  operation. 


Toble  B-1.  Elevation  differences   (mm)   relative  to  RANCH-D. 


Mmrk  1964-65 

1965-66 

1966-68 

1968-69 

1969-71 

1971-72 

A               -fl4 

+  1 

+  14 

-21 

+  6 

+  2 

B               +17 

+  5 

+  12 

-  9 

+  3 

+  2 

C               +10 

+  4 

+  7 

-12 

+  9 

-  1 

E                  -2 

+  7 

-  1 

_ 

-  4 

+  1 

F              +20 

-  1 

+  19 

-16 

-1 

+  7 

TEJON 

Figure  B-12a  shows  the  TEJON  (1964)  survey  figure.  There 
was  no  evidence  of  consistent  vertical  movement.  The  observed 
horizontal  direction  changes  between  the  three  surveys  suggest 
that  there  may  have  been  some  left-lateral  movement  at  a  rate 
of  1  to  2  mm/year  across  the  Oarlock  fault  (Parkin.  1966).  The 
taped  base  A-B  shortened,  introducing  some  compression  into 
the  figure. 

RANCH-TEJON  COMBINED  FIGURE 

Figure  B-1 3  shows  the  new  RANCH-TEJON  site.  The  level 
surveys  suggest  there  may  have  been  some  subsidence  of  marks 
K  and  NO  PORTAL  TUN  3  relative  to  B,  but  more  data  are 
needed  to  establish  trends.  The  results  of  the  honzontal  surveys 
indicate  a  general  contraction  of  the  figure  between  1970  and 
1972.  The  observed  direction  changes  hint  at  small  nght-lateral 
movement  across  the  North  Oarlock  fault  and  larger  left-lateral 
movement  across  the  South  Oarlock  fault.  The  position  vectors 
of  the  stations  north  of  PJW19  relative  to  PJW19  are  well  within 
their  error  ellipses,  while  the  vectors  for  RANCH  D  and  EMIL 
are  barely  within  the  error  ellipses.  These  results  suggest  a  con- 
tinuing rate  of  8  mm/year  movement  across  the  South  Oarlock 
fauh  (Miller,  1972d). 

It  should  be  pointed  out  that  measurements  on  line  51  of  the 
California  Geodimeter  Network,  which  crosses  the  Oarlock 
fault  near  RANCH-TEJON.  suggest  a  long-term  left-lateral 
movement  rate  of  5  ±  2  mm/year  (Savage  and  others,  1973,  p. 
48,  figure  3).  The  Geodimeter  results  are  consistent  with  the 
results  from  the  RANCH-TEJON  surveys.  This  indicates  that 
the  Oarlock  fault  in  this  area  may  be  more  active  than  previously 
suspected. 

QUAIL 

This  site  was  established  on  the  San  Andreas  fault  (figure 
B-14)  in  1964  and  resurveyed  in  1965.  1968,  and  1971.  The 
triangle  closures  in  1964  and  1965  were  somewhat  larger  than 
one  would  expect.  The  astronomic  azimuth  ofline  B-F  changed 
by  about  1.5  seconds  between  1964  to  1965,  and  1965  to  1968, 
but  the  value  in  1971  was  essentially  the  same  as  that  in  1964. 
The  taped  base  B-F  changed  by  4  mm  over  the  course  of  the 
surveys.  First-order  levels  were  run  in  1964,  1965,  and  1968. 

There  does  not  appear  to  be  any  consistent  pattern  to  the 
elevation  differences.  There  is  no  indication  of  fault  movement 
at  this  site  based  on  the  horizontal  survey  results  (Miller. 
1971b). 

HUGHES 

This  site  was  established  on  the  San  Andreas  fault  (figure 
B-1 5)  and  was  reobserved  in  1965  and  1971.  The  survey  results 
were  within  specifications,  and  the  astronomic  azimuths  and 
taped  distances  remained  essentially  constant.  First-order  levels 
were  run  in  1964,  1965.  and  1971. 

There  does  not  appear  to  be  any  consistent  pattern  in  the 
elevation  changes,  and  there  is  no  evidence  of  fault  movement  in 
the  vicinity  of  this  site  based  on  the  survey  results  (Miller, 
1972c). 

PALM 

This  site  was  established  on  the  San  Andreas  fault  in  1964  and 
resurveyed  in  1965  and  1966.  First-order  levels  were  run  at  the 
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site  in  1964  and  1965.  There  was  some  evidence  of  rotation  of 
the  astronomic  azimuth  of  the  base  line.  The  taped  length  of  the 
base  increased  by  8  mm  between  1964  and  1966.  The  tnangle 
closures  and  angle  corrections  were  quite  large  (Miller,  1966a). 

The  horizontal  surveys  revealed  no  displacement  across  the 
San  Andreas  fault  larger  than  the  survey  errors  (Miller,  1966a). 
The  level  survey  revealed  a  general  uplift  of  the  figure. 

WARM 

This  site  was  established  in  1964  on  the  Clearwater  fault  and 
resurveyed  in  1965.  The  1965  results  show  no  evidence  of  move- 
ment (Parkin,  1965b). 

CAST 

This  site  was  established  in  1964  and  resurveyed  in  1965,  1967, 
and  1970.  It  was  supposed  to  cross  the  San  Gabriel  fault,  but 
apparently  the  fault  was  missed  when  the  figure  was  laid  out 
(Miller.  1971a).  Therefore,  no  further  discussion  will  be  present- 
ed. 

BARREL 

This  site  on  the  San  Andreas  fault  was  established  in  1964  and 
reobserved  in  1965.  1966,  1967,  1970,  and  early  1971.  The  1970 
survey  was  made  between  December  1970  and  January  1971, 
and  the  1971  surveys  were  made  in  February  and  May  following 
the  1971  San  Fernando  earthquake.  The  astronomic  azimuth  of 
line  B-D  seems  to  be  3  to  6  seconds  too  large  in  1964.  The  taped 
length  of  B-D  varied  by  as  much  as  5  mm  between  the  surveys. 
Triangle  closures  and  corrections  to  angles  were  within  specifica- 
tions (Miller,  1971c). 

Figure  B- 1 6  shows  the  survey  figure  and  its  relation  to  the  San 
Andreas  fault  zone,  and  table  B-2  shows  the  elevation  changes 
relative  to  mark  B.  Note  that  the  total  net  change  in  elevation 
shows  that  all  the  marks  subsided  relative  to  B.  However, 
between  1965  and  1966,  the  marks  all  rose  relative  to  B.  The 
results  of  the  horizontal  surveys  show  that  no  fault  movement 
larger  than  the  survey  error  occurred  on  the  San  Andreas  fault 
near  this  site  (Miller  1971). 

Table  B-2.      Elevation  differences   (mm)   relative  to  BARREL-B. 


Table  B-3.      Elevation  differences   (mm)   relative  to  PEAR-G. 


Mark 

;p<M-dj 

1965-66 

1966-67 

1967-71 

A 

-4 

+3 

-2 

-1 

C 

-2 

-1-3 

-3 

-3 

D 

-1 

+s 

-2 

-2 

E 

-6 

+* 

-2 

-1 

F 

-10 

+  6 

-2 

-2 

PEAR 

This  survey  figure  was  established  across  some  minor  faults 
associated  the  San  Andreas  fault  zone  in  1964  (figure  B-17)  and 
reobserved  in  1965,  1966,  1967,  and  1969.  The  astronomic 
azimuth  of  line  A-G  changed  by  about  4  seconds  between  1964 
and  1965,  but  it  remained  essentially  constant  after  1965.  The 
taped  length  of  A-G  changed  by  6  mm  between  1964  and  1965. 
After  that,  changes  were  ±  2  to  4  mm  in  the  length  of  A-G.  The 
triangle  closures  and  angle  corrections  were  large  in  the  first  two 
surveys  and  much  smaller  in  the  last  two  surveys.  The  results  of 
the  horizontal  and  vertical  surveys  indicate  that  no  systematic 
movement  larger  than  the  survey  errors  has  occurred  (Miller, 
1969). 


Mark 

196*-65 

1965-66 

1966-67 

1967-69 

A 

-1-2 

0 

-2 

+2 

B 

-1-4 

-3 

-h3 

0 

C 

-1-3 

-1 

0 

0 

D 

-1-2 

-3 

-1-2 

-n 

F 

-1-2 

-3 

-1-1 

0 

H 

-1-1 

-1 

-1 

0 

WRIGHT 

This  site  was  established  on  the  San  Andreas  fault  in  1964  and 
resurveyed  in  1965.  The  survey  errors  were  somewhat  larger 
than  expected.  There  does  not  appear  to  be  any  systematic  move- 
ment on  the  fault  larger  than  the  survey  error  (Parkin,  1965a). 

CLEG 

This  figure  was  established  on  the  Cleghom  fault  in  1965  and 
resurveyed  in  1966  and  1967.  First-order  levels  were  run  in  1966 
and  1967.  Neither  the  horizontal  nor  the  vertical  surveys  indicat- 
ed any  significant,  systematic  fault  movement  (Anonymous, 
1967). 

CEDAR 

This  site  was  established  on  the  Cleghom  fault  in  1964  and 
resurveyed  in  1964  and  1967  (Miller,  1967b).  The  site  has  been 
destroyed  since  1967  (Greensfelder,  1973).  There  was  some 
movement  of  one  of  the  marks  in  the  figure.  Measurements  not 
involving  that  mark  indicated  no  fault  movement  (Miller, 
1967b). 

DEVIL 

This  site  was  established  on  the  San  Andreas  fault  in  1964 
(figure  B-18)  and  resurveyed  in  1965,  1966,  1968,  1970,  and 
1972.  Levels  were  run  in  1964,  1965,  1968,  1970,  and  1971.  There 
were  large  changes  (  ±  4  seconds  maximum)  in  the  astronomic 
azimuth  of  line  B-D  between  the  surveys;  the  taped  length  of 
B-D  has  increased  by  1  mm.  The  triangle  closures  and  angle 
corrections  were  generally  well  within  the  specifications  of  the 
survey.  The  results  of  the  horizontal  surveys  indicate  no  fault 
movement  between  1970  and  1972  (Miller,  1972b).  The  vertical 
surveys  indicate  that  the  northeast  side  of  the  San  Andreas  fault 
is  moving  up  relative  to  the  southwest  side.  The  vertical  surveys 
will  be  discussed  in  more  detail  under  "Survey  Results  Near  San 
Bernardino". 

RIALTO 

This  site  was  established  in  1964  (figure  B-19)  on  the  San 
Jacinto  fault  and  reobserved  in  1965.  1967.  1968.  1970.  and  1972. 
The  astronomic  azimuth  of  line  B-D  remained  essentially  con- 
stant until  the  survey  of  1972.  when  it  decreased  about  3  seconds. 
The  taped  length  of  B-D  varied  over  a  range  of  5  mm.  Triangle 
closures  and  angle  corrections  were  within  the  specifications  of 
the  surveys.  There  is  reason  lo  believe  that  mark  C  was  disturbed 
by  heavy  equipment  between  1968  and  1970  (Miller.  1972b). 

Levels  were  run  at  this  site  in  1964.  1965,  1967,  1968,  1970, 
1971,  and  1972.  A  level  line  connecting  RIALTO  to  COLT  was 
established  in  1971  and  reobserved  in  1972. 

The  results  of  the  horizontal  surveys  showed  a  right-lateral 
creep  rate  of  about  2  mm/year  up  until  1968.  Between  1968  and 
1970,  there  may  have  been  some  small  amount  of  fault  move- 
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ment;  but  the  rate  was  not  as  large  as  the  previous  rate.  From 
1970  to  1972,  there  was  no  fault  movement  larger  than  the 
survey  error.  It  is  interesting  to  note  that  the  elevation  changes 
also  changed  trend  in  the  period  1968  to  1970.  This  will  be 
discussed  in  more  detail  under  "Survey  Results  Near  San  Ber- 
nardino". 

COLT 

This  site  was  established  on  the  San  Jacinto  fault  in  1964 
(figure  B-20)  and  reobserved  in  1965,  1967,  1969,  1970,  and 
1972.  The  astronomic  azimuth  of  line  B-D  decreased  about  15 
seconds  from  1964  to  1972.  The  taped  length  B-D  increased  3 
mm  between  1964  and  1972.  There  is  some  indication  that  mark 
F  was  disturbed  between  1970  and  1972,  and  calculated  position 
vectors  suggest  that  marks  F  and  D  have  been  unstable  since 
1964  (Miller,  1972b).  Levels  were  run  at  this  site  m  1964,  1965, 
1967,  1968,  1969,  1971,  and  1972.  A  level  hne  connecting  this 
site  to  RIALTO  was  established  in  1971  and  reobserved  in  1972. 

The  results  of  the  horizontal  surveys  indicate  that,  between 
1964  and  1970,  there  was  fault  movement  with  a  rate  of  about 
1.5  mm/year  on  the  San  Jacinto  fault.  In  1972,  all  the  marks 
moved  southeast  relative  to  mark  A.  The  vertical  surveys  also 
indicate  a  substantial  change  in  1970.  These  results  will  be  dis- 
cussed in  more  detail  under  "Survey  Results  Near  San  Bernar- 
dino". 

SURVEY  RESULTS  NEAR  SAN  BERNARDINO 

In  this  section,  the  results  from  the  three  survey  figures  near 
San  Bernardino  will  be  discussed  in  more  detail,  and  an  attempt 
will  be  made  to  correlate  these  results  with  other  data.  First,  let 
us  recall  the  results  mentioned  previously. 

At  DEVIL,  there  was  no  indication  of  fault  movement  prior 
to  1970;  while  between  1970  and  1972,  there  was  15  mm  of 
right-lateral  movement  at  the  site.  At  RIALTO,  there  was  more 
or  less  steady  right-lateral  movement  of  2  mm/year  (total  about 
10  mm)  between  1964  and  1968  with,  perhaps,  a  slight  bit  of 
right-lateral  movement  between  1968  and  1970.  After  1970, 
there  was  essentially  no  movement.  At  COLT,  there  was  move- 
ment with  a  rate  of  about  1.5  mm/year  (total  about  9  to  10  mm) 


between  1964  and  1970,  and  there  was  no  evidence  of  movement 
after  1970.  The  level  data  at  DEVIL  show  a  more  or  less  steady 
uplift  of  the  northeastern  block  relative  to  the  southeastern  block 
relative  to  the  southwestern  block,  but  the  data  are  too  sparse  to 
draw  firm  conclusions  about  possible  changes  in  rates.  At  RI- 
ALTO, the  level  data  indicate  a  change  in  trend  between  the 
1968  and  1970  surveys  (figure  B-20a);  while  at  COLT,  a  change 
in  trend  occurs  between  the  1970  and  1972  surveys.  Thus,  both 
horizontal  and  vertical  data  at  COLT  and  DEVIL  suggest  a 
major  change  in  the  strain  field  between  1970  and  1972. 

Three  Geodimeter  lines  of  the  California  Geodimeter  Net- 
work have  been  measured  near  this  area;  they  are  lines  66,  67, 
and  70  (Greensfelder  and  Bennett,  1973,  figure  6).  Lines  66  and 
67  show  length  increases  of  4  to  5  cm  between  1968  and  1970 
and  length  decreases  of  about  2  cm  between  1970  and  1972 
(figure  B-21).  Line  70  lengthens  by  about  10  cm  between  1968 
and  1972.  Thus,  the  Geodimeter  data  also  suggest  a  change  in 
the  strain  field  between  1970  and  1972.  If  the  hypothesis  of 
periodic  fluctuations  in  line  lengths  proposed  by  Greensfelder 
and  Bennett  (1973)  is  vahd,  the  data  from  these  small  survey 
figures  may  be  reflecting  this  change.  More  data  are  needed  to 
resolve  this  problem. 

It  is  interesting  to  note  that  on  September  12,  1970,  a  magni- 
tude 5.4  earthquake  occurred  northwest  of  the  San  Bernardino 
area  and  that  on  February  9,  1971,  the  magnitude  6.6  San  Fer- 
nando earthquake  occurred  100  km  east  of  the  San  Bernardino 
area  (figures  B-19b,  B-20b,  and  B-21).  It  is  possible  that  the 
change  in  the  strain  field  between  1970  and  1972  was  related  to 
one  or  both  of  these  two  events,  but  the  data  are  not  adequate 
to  resolve  this  question. 

Finally,  it  is  intriguing  to  note  that  Kanamori  and  Chung 
(1973)  and  Wu  and  others  (1973)  found  changes  in  seismic 
velocities  in  the  vicinity  of  Riverside  between  1966  and  1971. 
They  indicate  that  these  changes  might  be  related  to  future 
earthquake  activity  in  the  vicinity  of  Riverside.  This  leads  us  to 
the  possibility  that  the  small  survey  figures  may  be  giving  us 
information  on  dilatancy  preceding  an  earthquake.  More  fre- 
quent measurements  on  somewhat  expanded  figures  and  level 
nets  should  be  made  to  test  this  important  possibility. 
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Figure  B-2.        Plan  of  the  survey  figure  at  the  MAR  site. 


Figure  B-3.        Plan  of  the  survey  figure  at  the  STADIUM  site. 


PFS 


Figure  B— 4.        Plan  of  the  survey  figure  at  the  CAMP  PARKS 
site. 
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Figure  ^5.        Plan  of  the  lurvey  figure  at  the  UNION  site. 


Figure  B— 6.        Plan  of  ttie  survey  figure  at  the  VERAS  site. 
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Figure  B-7.        Plan  of  the  survey  figure  at  the  OS- 
GOOD site. 
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Figure  ft-8.       Plan  of  the  survey  figure  at  the  TEM  lite. 


Figure  ft-9.       Plon  of  the  survey  figure  at  the  SANTA  site. 
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Figure  B-10.        Plan  of  the  survey  figure  at  the  MEADE  site. 
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Figure  B-1 1.        Plan  of  the  survey  figure  at  the  RANCH  site. 


Figure  &-12o.   Plan  of  the  survey  figure  at  the  original 
TEJON  site  (1964-66). 
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Figure  B-1 2b.   Plan  of  the  survey  figure  ol  the  new 
TEJON  site   (1V70). 
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Figure  B-13.      Plan  of  the  survey  figure  at  tfie  RANCH 
-TEJON  site  (1970). 


Figure  ^14.      Plan  of  ttie  survey  figure  at  the  QUAIL 
site. 
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Figure    B-15.      Plon    of    tf>e    survey    figure    at    the 
HUGHES  site. 
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Figure  &-21.      Geodinwter  line  length  changes,  Polmdale  to  San  Bernardino  (from  Greenifelder  and  Bennett,  1973). 
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